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1. INTRODUCTION 


Our ability to undarstand and modal turbulent flows still 
relies heavily on the availability of accurate measurements 
of mean and fluctuating quantities within the flow. Until recently, 
the hot wire was the only reliable tool available for the measurement 
of fluctuating velocities. In fact, almost all of our present 
knowledge about turbulent flows is based on measurements made 
with hot wires. In relatively simple flows (moderately two- 
dimensional with small cross-flows), reliable and accurate hot 
wire measurements are now possible with fully automated data 
acquisition and reduction systems which minimize errors due 
to drifts in calibrations. An example of such a system is given 
in Ref. 2. 

However, as we turn our attention towards more complex 
turbulent flows, a need for more sophisticated measurement techniques 
has become apparent. These complex flows include those with 
compressibility affects, strong three-dimensionality (with steep 
mean gradients), flow reversals, and time-dependent behavior. 
Since about the mid-sixties, the most popular alternative tool 
for measuring mean and fluctuating velocities in turbulent flows 
has been the Laser Doppler Velocimeter (LDV) . 

The most popular LDV arrangement used for wind tunnel measure- 
ments is the dual beam or fringe method. In this method, one 
of the laser lines is split into two lines of equal intensity 
which are then focused through a lens so that thay cross over 
at the focal point. The flow is seeded with small particles 
(typically less than 3-4 )m in diameter) which follow the fluid 
motion. A3 these particles pass through interferometric fringes 
created by the crossed laser beams, light is scattered off them 
which is received by a photodetector . The frequency of this 
scattered light, along with a knowledge of the fringe pattern 
formed by the lasor beams, provide the means to calculate the 
velocity of the partible. The fringe method, especially in 
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the forward scatter, off-axis mode, generally offers the best 
signal-to-noise ratios and spatial resolution. 

Although LDV systems are somewhat complex and tiresome 
to set up, they have certain advantages over hot wires for turbulence 
measurements. The fact that Laser Doppler Velocimetry is non- 
intrusive is especially beneficial in the measurement of unstable 
flow phenomena which are very sensitive to the presence of measure- 
ment probes. In certain situations, LDV systems can also provide 
greater spatial resolution and better directional descrimination 
than hot wires. This makes it possible to use LDV systems for 
the measurement of separated flows. Since an LDV measures the 
velocity directly, independent of the thermodynamic properties 
of the flow, it is particularly attractive for velocity measurements 
in compressible flows. Furthermore, the calibration converting 
the frequency to velocity is linear and easy to implement in 
software. This feature also allows for uniform sensitivity 
in measuring both moderate and high turbulence intensities. 

Two-color LDV systems capable of measuring two components 
of velocity simultaneously are now being widely used. However, 
the main interest in the present investigation was to study 
three-dimensional interactions where it is desirable to obtain 
measurements of all three velocities. Hence, the first objective 
was to develop a laser velocimeter system capable of measuring 
all three components of velocity simultaneously so that all 
six components of Reynolds stress may be computed. Another 
objective was to compare these LDV measurements directly with 
those obtained with hot wires in flow fields where both techniqes 
are expected to perform satisfactorily. The purpose of this 
is to evaluate the performance of the system quantitatively 
and objectively. 

The approach followed in the present investigation was 
to convert an existing two-component LDV system into one capable 
of measuring all three components of velocity simultaneously. 
The system utilizes two wavelengths (488.0 and 514.5 nm) from 
a 4 -watt Argon-Ion laser. The main four-beam matrix measures 
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u and v directly. The green line in the four-beam matrix is 
split (in half) using a dichroic filter and directed over the 
top of the traverse mechanism with mirrors, giving the third 
beam pair for the measurement of the w-component. This pair 
of beams (with rotated polarization) measures w sin 45° + v cos 4^. 
Scattered light is collected in the off-axis forward scatter 
mode using two collection lenses. 

Signal processing is accomplished with single-particle 
burst counters, and the validated data are multiplexed through 
a "home-built" interface to an HP 9845B desk-top computer. 
Some selected first and second order products are reduced on-line, 
and the raw data is dumped onto floppy disk. An off-line program 
reduces the data, giving up to third order quantities and also 
plots histograms of the raw data for each channel. The software 
includes the capability to filter out noise by examining the 
histograms . 

The optical and signal processing hardware is described 
in Section 2. The data acquisition and reduction software is 
described in Section J and detailed operating procedures are 
given in Section 4. Some problems, inherent to 3-component 
LDV systems, are preserved in Section 5. Sample results from 
an experiment measuring mean and turbulence quantities in a 
vortex/mixing layer interaction are compared directly to results 
obtained using crossed hot-wire anemometry in Section 6, and 
concluding remarks are presented in the fina~ section. Complete 
software listings written in BASIC and ASSEMBLY languages to 
run on the HP 9845B desk-top computer are included in the appendix. 

2. OPTICAL SYSTEM AMD SIGNAL PROCESSING HARDWARE 

The hardware for the 3 -component LDV system can be divided 
into three categories: the optical system, the signal processing 
instrumentation, and the computer. The LDV optics consist of 
the optics table, where the laser beam is split into green and 
blue beam pairs, the transmitting optics, where the beams are 
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directed into the flow field, and the receiving optics, where 
the scattered light is picked up by photodetectors. The signal 
processing instrumentation consists of nplifiers, filters, 
burst counters and a computer interface. ’ . sampling procedures 
are all computer controlled. 

2 . 1 Optics Table 

The optics table consists of the laser and all the optical 
elements needed to provide the necessary four-beam matrix. 
Fig. 1 shows a schematic of the arrangement with component numbers 
as referred to in this section. A 4-watt Argon-Ion laser (Lexel 
Model 95) is used to produce the main beam. The beam then passes 
through a collimator (1), which ensures that the beam waist 
occurs at the focal point of the transmitting lens. 

The collimated beam is passed through a color separator 
box, which consists of a polarization rotator (2) , an attenuator 
(3), and a pair of high dispersion Brewster angle prisms (4) 
which are used to separate the multi-line beam into two colors, 
blue (488 nm) and green (514.5 nm) . These two beams are then 
reflected across the box by mirror (5) , and out of the box by 
mirrors (6, 7). 

Follow ing t hg color separation, the green beam's polarity 
is rotated to horizontal (8) , and the beam is split into two 
beams in the vertical plane (9) . Most beam splitters prefer 
this type of perpendicular polarization for maximum efficiency. 
Using the beam displacer (10), the blue beam is then moved to 
the center of the optics, and its polarity is rotated to vertical 
(11). The blue beams are split in the horizontal plane (12). 
At this point, the four beams are each displaced 25 mm from 
the optical axis. 

Two Bragg cells (13, 14) are used to shift the frequency 
of one beam from each pair. The unshifted beam passes through 
an optical rod so that the path lengths are matched. The frequency 
is shifted by a fixed amount of 40 Mhz. This shift creates 
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a moving systam of fringes at the beam intersection point, allowing 
for directional discrimination of the velocity. Frequency shifting 
also helps to reduce the percentage of frequency change in highly 
turbulent flow, to reduce fringe bias, and to optimize frequency, 
thus enabling easy removal of the pedestal by high pass filtering. 

The four beams are then passed through a beam steering 
module (15) . The module consists of a set of wedge prisms that 
can be independently rotated about the beam axis to steer the 
shifted beam in any direction. This allows a more precise alignment 
of the beams. The beams are finally passed through a beam displacer 

(16) to reduce the beam spacing to 13 mm and a rotating prism 

(17) before leavinq the optics table. The rotating prism enables 
the four-beam matrix to be rotated independently so that the 
beams may be aligned relative to the tunnel axes. 

2 . 2 Transmitting Optics 

The transmitting optics (Fig. 2) are mounted on a traversing 
mechanism with three degrees of freedom. The traverses are 
driven by individual stepper motors. The four beams from the 
optics table are directed by a set of five mirrors through a 
dichroic filter before being focused by a 380 mm (15 inch) focal 
length lens. This main four-beam matrix measures u and v directly. 
The dichroic filter, set at an angle of about 30 degrees to 
the incoming beams, splits the green beams in half, which provides 
the third beam pair for the third velocity component, w. Thi3 
third beam pair is directed over the top by mirrors, passed 
through a polarization rotator (giving it a different polarity 
than the main-axis green beam pair) and then focused at the 
focal point of the main beam set. Since this third beam pair 
intersects the main measuring volume at a 45° angle to the main 
axis, it measures v and w with equal sensitivity, with the measured 
component being w v ■ (v cos 45° + w sin 45°) . Since v i3 measured 
directly, w can be evaluated using the equation: 
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w. , - v cog 45° 

w - sin 45° (1) 

Typical p-obe volume dimensions for each beam pair in the 
present configuration are 10 mm in length and 0.2 mm in diameter 
(Fig. 3) . However, the actual "viewed" dimensions are reduced 
considerably, as discussed below, in Section 2.3. 

2.3 Receiving Optics 

The detector system (Fig. 4) is in the off-axis forward 
scatter mode. The receiving optics are mounted on a traversing 
gear, also run by stepper motors, which moves synchronously 
with the transmitting optics. Scattered light is collected 
by two 380 mm (15 inch) focal length lenses. The collimated 
light is passed through filters to separate the colors and the 
off-axis line is additionally passed through a polarization 
filter to avoid collecting light scattered by the main axis 
green pair. The collected light is then focused by 250 mm (10 
inch) focal length lenses onto pin-hole apertures mounted in 
front of the three photomultipler tubas. The collection angle 
and diameter can be adjusted to select the effective (viewed) 
probe length (Fig. 5) . In the present set-up, a collection 
angle of 30 degrees and an aperture diameter of 0.5 mm were 
used to give an effective length of about 1.5 mm. 

2.4 Signal Processing Hardware 

The signals from the photomultiplier tubes are amplified 
and relayed to the signal processors via high-pass filters and 
mixers (Fig. 6) . The amplifiers used are EIN model 403LA with 
a fixed gain of 37 dB and the filters are Allen Avionics F2440 
with a fixed high-pass cutoff of 10 MHz. The filtered signals 
are mixed electronically with sine waves from three Tektronix 
SG503 levelled Sine Wave Generators. The mixer is commercially 
available from Hewlett-Packard, model 10534A. The mixing procedure 
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is necessary for low-speed flows, where the actual Doppler fre- 
quencies are small compared to the Bragg frequency of 40 MHz. 
So in order to reduce the effective measured frequency and hence 
improve the counter resolution, the incoming signals are mixed 
with sine waves of known frequency. 

The mixed signal is fed into single particle burst counters 
(Macrodyne model 2C96-2 and 3003) via high-pass/low-pass filters 
and an amplifier (xlO). The counters measure the zero crossings 
of the Doppler signal, which is related to the Doppler frequency 
by the range set on the Macrodyne. The range is set manually 
based on expected flow velocities, since it limits the frequencies 
that the processor can see for the given 10 bits of resolution. 

The processors use two checks to validate a Doppler signal. 
The first check is the usual 5/8 comparison, where the processor 
checks the frequency for 5 zero crossings against that for 8 
crossings. The second check is the multi-sequence check. Positive 
and negative thresholds are set on the signal, and a validated 
output is permitted only if, for all eight fringe crossings, 
the signal passes through a positive threshold, a zero level, 
and a negative threshold in the proper sequence. The digital 
data (consisting of a 10 bit data word with 3 bits giving the 
range) and a sync pulse (produced every time the front end of 
a valid burst is detected) are passed to the computer interface. 

2 . 5 Computer Interface 

A NASA LDV-A/D computer interface (Cl) is used to transfer 
data from the LDV signal processor to the computer (Fig. 7). 
The Cl can interface either digital or mixed analog and digital 
data to an HP 9845B desk-top computer. The Cl consists of an 
eight-channel multiplexer, a four-channel A/D converter, and 
an event synchronizer with time interval counter. 

For use with the LDV, tne inputs to the Cl are all digital. 
Six of the eight vord3 come from the processor, and two are 
time and status words from the synchronizer. The inputs are 
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multiplexed to a single digital data channel output. 

The Cl can accept data from the three processors in either 
random mode or sync mode. In random mode, the Cl will accept 
data inputs when an event occurs on any of the thrae channels. 
A dead time (between 5 and 50 jis) is set in this mode, which 
controls the minimum time between samples to ensure that a given 
particle is sampled only once. In sync mode, the Cl will accept 
data inputs when all three processors sample simultaneous events. 
In this mode, a coincidence time (between 5 and 50 ps) has to 
be set, which determines the time window within which all three 
events must occur. A detailed description of the Cl can be 
found in Ref. 5. 

• r * 

2.6 HP 9845B Desk-top Computer 

Multiplexed data are passed to the HP computer from the 
NASA >1DV-A/D using the HP 98032A high-speed 16-bit parallel 
interface. Jumpers labeled N S,B,D" are connected inside the 
98032A for proper operation with the LDV-A/D. A select code 
of 10 (screw setting on the 98032A) is set for use with the 
software described in the appendix. A data buffer of 24 kbytes 
is provided in the memory of the HP 9845B desk-top computer 
for storage of up to 2000 samples obtained from one measurement 
location. Since three of the six data words passed for each 
sample are merely monitored and discarded as described above, 
12 kbytes of raw data remain to be stored for each point. An 
HP 9895A floppy-disk drive is used for archival storage of raw 
data. The buffered data are written in real time to a sequential- 
access floppy-disk file. Enough header information is written 
to each file to identify the run, and to reproduce calibration 
tables. 

3. DAI A ACQUISITION AND REDUCTION SOFTWARE 

Data aquisition and reduction on the HP 984 5B is dona via 
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two programs. Tha data acquisition program controls the processor 
sampling. The program accepts heading arid initialization parameters 
provided by the operator, performs a fast I/O handshake to acquire 
the raw data from the LDV-A/D computer interface, and writes 
the raw data and initialization parametars onto floppy-disk. 
In addition, the data acquisicion program has a limited capability 
for on-line data reduction so that key results may be monitored 
during a run. The off-lir.e data reduction program converts 
the raw data into instantaneous velocities and computes al’ 
the statistical quantities. The off-line program also has options 
to plot histograms of tha raw data, to plot profiles of the 
reduced data, and to filter out noise. 

3 . 1 Data Acquisi tion Program 

The data acquisitic program structure can bo divided into 
four areas: (1) initialization of variables, (2} data acquisition, 
(3) lini-ed data reduction, and (4) raw data storage ,sse tha 
block diagram in Fig. 8) . During the initialization stage of 
the program, relevant- test parameters are input to be saved 
along with the data. Tr-asa parameters include the spatial coord- 
inates, Bragg and mixing f requencies , and angle of the third 
beam pair. Fringe spacings and probe dimensions are also calcu- 
lated at this time. 

Data i3 acquiiad from the LDV computer interface by mean? 
of a fast I/O handshake. The operator requests tnat a certain 
number of data samples be taken and stored in the data Luffei 
as described in Section 2.6. The 16-bit data word3 are put 
into an integer array, three words wide, containing values for 
wha three velocity components The number acquired from the 
counter occupies only bits 0-2 of the data word. Bits 10-13 
specify the range, and bits 14 and 15 are unused. The data 
array, along with the initialization parameters given by the 
operator, are stored in a floppy-disk file. An an ontior some 
of tha data can be reduced on-line, before tne raw data is written 
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to floppy-disk. This on-line reduction feature uses the same 
algorithms as the off-line data reduction program, but is faster 
since it computes fewer guantitites. Data acquisition timns 
depend on many factors other than the software, but on-line 
data reduction takes about 30 seconds per point on the HP 9845B, 
and each file write to the floppy-disk requires about 20 seconds, 
assuming 2000 samples are taken per point. 

3.2 Data Reduction Program 

A separate off-line program is used to reduce the complete 
data. The data reduction program has four major sections: 
(1) data read from floppy-disk, (2) conversion of raw data to 
instantaneous velocities, (3) calculation of sum3 of instantaneous 
velocity components, and (4) calculation of moments from the 
sums (see the block diagram in Fig. 9) . Since the HP 9845B 
is a relatively slow micro-computer, the bulk of the calculation 
routines are written in ASSEMBLY language to reduce computation 
time. 

The program reads the initialization parameters and raw 
data which were written on floppy-disk by the data acquisition 
program. The raw data is used together with the mixing frequencies, 
Bragg shift frequency, third beam crossing angle, and fringe 
spacings to calculate the three velocity components for each 
sample. The velocities for each sample are calculated in turn, 
and a running sum is maintained. After the samples have been 
summed, the average values of the various moments are computed, 
and the mean velocities and turbulence quantities are calculated 
from them. Reading from the floppy-disk takes about 20 seconds, 
while the data reduction including histogram plots requires 
about 90 seconds per point on the HP 9845B. 

The procedure for data reduction is fairly straightforward. 
The raw data is converted to frequencies using the range 3et 
on the processor. The measured frequency is given by the relation 
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f rea - 3.2 * 10 4 / 2 R * D 

where D - raw data. 

The velocities are then calculated from the resultant frequen- 
cies, by the relation: 

velocity - d f * ((ffaragg- ^mix)” «res) (3) 

Once the velocities have been evaluated, average values of the 
various moments are computed as defined below: 


Su - S<u(i)> / Ns 
Sv - S<v(i) > / Ns 
Sw « S<w(i)> / Ns 
Suu - S<u(i)u(i)> / Ns 
Suv ■ S<u(i)v(i)> / Ns 
Suw « S<u(i)w(i)> / Ns 
Sw - S<v(i)v(i)> / Ns 
Svw ■ S<v(i)w(i)> / Ns 
Sww ■ S<w(i)w(i)> / Ns 
Suuv » S<u(i)u(i) v(i) > / Ns 
Suuw - S<u(i)u(i)w(i) > / Ns 
Suw - S<u(i)v(i}v(i)> / Ns 
Suww » S<u(i)w(i)w(i) > / Ns 
Suvw - S<ufi)v(i)w(i)> / Ns 


( 4 ) 


Using these definitions, the signal statistics are then calculated 
assuming nearly infinite sample size: 

u - Su 

v « Sv 

w - Sw 
2 

U 1 - SUU - SuSu 
2 

v ' « Sw - SvSv (5) 
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The implementation of the data reduction software is somewhat 
more complex. In order to achieve reasonable running times, 
the ASSEMBLY code is optimized for speed rather than clarity 
of operation. Operations which do not change between samples 
are performed only once. Thus, 

veJocity - d f * ((fbragg ~ f mix) “ f res) * 1q€ («) 

becomes 

f int " *bragg “ f mix ; d int " d f * 106 ( 7 ) 

hence, velocity - d int * (fi n t “ f res) (8) 

Additional calculations are required to obtain the w component 
velocity. Since the third set of beams measures w v - w sin 45° 
+ v cos 45° , to find w, we must also perform the calculation: 
w » (w v - v cos 45°) / sin 45°. (9) 

The raw data output of the digitizer is in the form of 
a 10 bit data word, giving a range of possible values (counts) 
from 0 to 1023. Since 2000 samples are normally taken for each 
point, a given value may be encountered many times. Accordingly, 
each time a new value of the raw data is encountered, the cor- 
responding velocity is calculated and stored in a look-up table. 
The next tine that value is encountered by the program, the 
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proper velocity can be easily looked up, eliminating the need 
for another time-consuming real variable calculation. The exact 
running procedure of the data reduction program therefore, is 
as follows: 

1. Reads header and all raw data for a particular point from 
floppy disk. 

2 . Strips bits 10-13 from three .-aw data words corresponding 
to the three different velocity components, and uses them to 
calculate the three ranges set on the A/D. 

3. Calculates various intermediate values which ramain the 
same throughout the point. 

4. Reads a data word, strips off bits 0-9, and uses this raw 
datum as an index to look up it3 corresponding velocity. 

5. If it is a new value, the program calculates the velocity 
using the raw datum, and stores it in the appropriate place 
in the look-up table. 

6. For the w velocity component, it finds the actual velocity 
from w « (w v - v cos 45°) / sin 45 a . 

7. Performs steps 4, 5, and 6 three times— once for each velocity 
component of the sample. 

8. Updates the running sums of the velocity components and 
products of the velocity components. 

9. Performs steps 4 through 8 for all samples. 

10. Uses the sums to obtain the average velocities, Reynolds 
stresses, and third order products, and then prints these quantities. 
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11 plots histograms of the raw data for all three channels. 

12. Performs steps 1 through 11 for each profile point. 

13. Plots profiles of the reduced data. 

14. Tabulates normalized data profiles. 

15. Writes a summary file containing the reduced data to disk. 

The data reduction program also has a routine to filter 
noise ai\d spurious data. Each of the 2000 data samples has 
3 counts corresponding to the 3 velocities associated with it. 
The filtering routine causes the data reduction program to ignore 
samples associated with counts which are excessively far from 
the mean. 

The filtering routine first sorts the data into three frequency 
tables, one for each channel. In a table, each count, i, has 
associated with it a number, S^, which is the number of samples 
with that particular count. 

The filtering routine finds the average count for each 
channel by going through the frequency tables and using the 
formula 


where j - 1,2,3 ■ channel number 
Cj ■ average count 
Ns - number of samples. 

Next, the filtering routine finds the standard deviation by 
performing the summation 

1024 

CTj 2 - (11) 

Ns 
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where (7 j ■ standard deviation. 

An input variable called S^ QV is read by the filtering 
routine. All counts further than S<i QV standard deviations away 
from the mean will be filtered out. This is done by multiplying 
s dev b y <7j for each channel, and going through the frequency 
tables one more time. If the magnitude of (i - Cj) > Sa ev * (Jj, 
then Sjj is sat equal to zero. 

The routine which calculates the statistical quantities 
takes each sample one at a time. The three counts associated 
with each sample are found and looked up in the frequency table. 
If for any count the table entry (S^j) is zero, then the sample 
is discarded. Thus any sample which is excessively far from 
the mean in any one -of its three counts is not used. 

4 . OPERATING PROCEDURE 

4.1 Alignment Procedure 

To achieve the best possible beam crossing and the most 
effective measuring volume, each module in the optics system 
must be carefully aligned. Detailed alignment instructions 
for individual optical components are given in Ref. 3. The 
overall alignment procedure is described in this report. Component 
numbers in this section refer to those shown in Fig. 1. 

The first step is to check that the laser output beam is 
parallel to the optics table at the specified height, using 
the system alignment blocks. The collimator (1) should be positioned 
so that the laser beam goes through the center of the len3, 
and focused so that the beam waist occurs at the cross-over 
point. The aligned beam then passes through the polarization 
rotator (2) and into the color separator box. The attenuator 
(3) ensures that the beam has horizontal polarity at this point. 

Each component of the color separator must be aligned separ- 
ately. Tho dispersion prism (4) should be aligned so that the 
path length is the same in both prisms and is parallel to the 
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base of each \.»*lsm. Several beams emerge from the prism, with 
the two brightest ones being green (514.5 nm) and blue (488 
nm) . The two beams are reflected out of the box using mirrors. 
The beams should be centered on mirrors (6) and (7) . As the 
beams come out of the box, alignment blocks are used to check 
the beam positions. If both beams are off-axis in the same 
direction, mirror (5) is used to align them. If only one beam 
is off, the appropriate mirror, (6) or (7) , is used for the 
adjustment. Beam splitter efficiency is maximized when the 
beam polarization is perpendicular to the plane of the split 
beams. This is achieved through components (8) and (11). The 
beam splitters (9) and (12) cannot split the beam intensity 
exactly in half, so one beam of each pair is always slightly 
brighter. Since the Bragg cells (13) and (14) normally attenuate 
the shifted beam, the brightest beam is frequency shifted so 
that the output beam pairs have nearly equal intensities. Detailed 
instructions for aligning the Bragg cells are given in Ref. 4. 
After passing through the Bragg cells, directional wedges are 
used to project the beams onto a distant surface (~ 3 m) . Any 
beam misalignments are more easily seen this way, and with the 
use of a marked mask, the beams can be adjusted to the correct 
orientation. Mirrors (6, 7) are used to adjust the unshifted 
beams, and the beam steering modules (15) are used for rotating 
the shifted beams along two circular arcs. 

The next check is to ensure that the beams are parallel 
to each leg of the traverse mechanism on the transmitting optics 
table, so that beam alignment is maintained while traversing. 
The dichroic filter (Fig. 2) is adjusted so that the green beam 
pairs are split equally. Once the beams pass through the trans- 
mitting optics, the four beams must be arranged so that they 
all cross at the same position. A microscope objective is used 
to view the beam crossing. If the four beams are not symmetric 
about the optical axis, mirrors (6) or (7) can be adjusted to 
correct this. The beam steering modules (15) are used to ensure 
that the beams cross at the same point along the axis. The 
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third pair of beams are now aligned (by eye) so that they also 
cross at the sane point as the nain line beans. 

To align the receiving optics, a piece of translucent tape 
is placed at the bean intersection point to scatter the laser 
light. By tracking the scattered light, the receiving optics 
are aligned to give a sharp inage at the pin-hole aperature 
in front of the photonultiplier tubes. 


4.2 Signal Processing 
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The sensitivity of the photonultiplier (PM) tubes used 
in the receiving optics can be varied by varying the appliod 
voltage. Typically, a voltage of about 1000 volts is applied. 
This voltage can be increased to make the PM tubes nore sensitive, 
as long as the threshold levels on the processors are increased 
accordingly, since the anount of noise picked up is also in- 
creased. 

The neasured signal is nixed electronically with the signal 
from a sine wave generator. The frequency of the sine wave 
(the nixing frequency) is chosen based on the expected flow 
velocities. The nixing frequency is chosen such that the difference 
between it and the Bragg shift is about twice the naxinun expected 
Doppler frequency. This allows enough nargin for fluctuations 
about the expected Doppler frequency and still have a remaining 
nonzero resultant frequency (f res - *bragg “ f nix “ f dopp) ► 
If faix is too high, a biasing results, similar to the fringe 
biasing caused by stationary fringes. The nunber of fringes 
crossed by a particle per second (as seen by the processor) 
is proportional to At(fb ra gg - f a i x ) , where At is the time taken 
by a particle to cross one fringe. Noting that At is only determined 
by the fringe spacing and the flow speed, if f m i x i3 increased # 
the number of fringes crossed by a particle is effectively reduced. 
This means that signals from particles which cross the fringes 
at an angle may not have enough fringe crossings (8) to be validated, 
and hence a bias towards particles movirg perpendicularly to 
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the fringes (higher velocity) results. This gives a higher 
mean velocity but a lover fluctuation level. 

The amplitude of the sine wave must be chosen so that an 
adequate signal-to-noise ratio is maintained. Typically, a 
peak-to-peak amplitude of 1 volt is required. 

The mixed signal is then fed into single particle burst 
(Macrodyne) processors. The high pass and low pass filter frequen- 
cies are set so that the processor frequency is centered between 
the two. Usually the high pass filter is set at 0.5 Mhz, which 
is the lowest non-zero setting, and the low pass filter is set 
anywhere from 2 to 16 Mhz, depending on the magnitude of the 
velocity component being measured. The filter bandwidth should 
bo broad enough so that no parts of the fluctuating signal are 
attenuated. 

The gain on the processor is normally set to 10. The output 
of the processors is displayed on an oscilloscope, and the Doppler 
signals should read about 1 volt peak-to-peak. The PM tube 
voltage can be adjusted so that the signal is at the desired 
level. Signal levels of more than about 1 volt end up being 
clipped and will therefore not be validated by the processor. 

The comparator accuracy for the 5/8 signal validation test 
can be set between 0 and 10 count variation. The 0 setting 
is the most accurate, and 11 is off. The processor manufacturer 
recommends that this level be usually set to 9. 

The range on the processor sets the bandwidth of frequencies 
that the processor can see, according to Table 1. For each 
range, a particular frequency corresponds to a count, from C 
to 1023. The processor frequency should be matched to a number 
in the central column in Table 1, and the corresponding range 
should be set on the Macrodyne processor. An additional check 
can be made on the range setting by monitoring the analog output 
from the Macrodyne on a DC voltmeter. The output ranges from 
0 to 10 volts, which corresponds to the counts from 0 to 1023. 
The correct range setting is that which gives about 5 volts 
on the voltmeter at the operating velocity. 
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Tha threshold Is then set for the multi-sequence check. 
The threshold should be set so that the data rate is about half 
of the data rate at zero threshold. Another check is to block 
one beam in each color pair. Tha data rata on the corresponding 
channel should be zero when one beam is blocked. 

To obtain good data rates in air, the flow must be seeded 
with uniformly sized particles. Smoke, obtained from burning 
mineral oil or Incense, provides particles of approximately 
2 pm. 

The computer interface should be set up as shown in Fig. 10. 
First choose the coincidence mode. If shear stresses are to 
be evaluated, coincidence on all 3 channels must be selected. 
The coincidence time (5 to 50 ps) sets the window width during 
which coincidence is defined. This should be set at 5 ps. 
If only the individual velocities are desired, then the random 
mode can be selected. The dead time (5 to 50 ps) should then 
be set so that data from one particle is not recorded twice. 
A setting of about 25 ps is recommended for low-speed flows. 

The number of word3 that must be multiplexed can be calculated 
as follows: fwords - #inputs + 2. This n”nber should be rounded 
to the nearest even number. In the present case a setting of 
six is used. The event mode is set to LDV (digital data only) . 
The LDV- A/D switch en ables both dig ital and analog data to be 
interfaced simultaneously. The counter clock frequency is set 
to equal the approximate data rate, and the computer select 
is set to HP. 

5 . INHERENT PROBLEMS 

Some design problems inherent to 3-component LDV systems 
are discussed in this section. One main problem with some of 
the earlier designs which called for splitting a color to create 
the third beam pair had to do with cross-talk. This is w.iere 
signals from tha two channels bearing the same color could not 
be adequately separated. In the original design of the present 
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set-up, the polarization of one of the green beam pairs was 
rotated relative to the other, so that the receiving optics 
could distinguish between the two signals. However, using a 
relatively large angle for tne off-axis beams (4 5° ) and two 
separate collection lenses, cross-talk between the two green 
channels has been almost eliminated, thus making the polarization 
rotation dispensable. 

Some earlier designs of 3-component LDV systems measured 
the u + w velocity component with the off-axis third beam pair. 
It is shown in Ref. 1 that measuring the v + w component instead, 
as done in the present system, reduces the uncertainty in the 
w component relative to the uncertainty in w from these earlier 
systems . 

Another problem has to do with signal coincidence. Details 
of the probe volumes for the present system are shown in Fig. 3. 
It is clearly illustrated how the cross-over region between 
the three sets of beams forms a very small fraction of the overall 
probe volume. Thus, with heavy seeding (necessary for three-channel 
work) , the electronics may validate data received from different 
particles which are not necessarily in the cross-over region 
but are within the coincidence time window set on the interface. 
This results in a lack of correlation between the measured veloc- 
ities, and causes the evaluated shear stresses to be inaccurate. 
Two schemes have been used in an attempt to minimize this problem. 
First, the coincidence time was made so short that measurements 
from different particles may be considered instantaneous. The 
minimum setting of 5 jis available in the present hardware was 
used; this is equivalent to about half the flight time of a 
particle passing through the probe volume. The second procedure 
involved reducing the effective "viewed” probe length and thereby 
reducing the probability of thi3 "apparent coincidence" (Fig. 5) . 
A collection angle 0 of 30° and an aperture diameter of 0.5 mm 
were used to give an effective length of about 1.5 mm. 
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As an initial check on the accuracy of the present system, 
some preliminary measurements have been made in a vortex/mixing 
layer interaction, previously investigated using the cross-wire 
technique (Ref. 6). Since the induced cross-flow angles in 
this interaction are only 5° -10° , cross-wire measurements are 
expected to be accurate to within about 5%. A schematic of 
the experimental set-up is shown in Figs. 11a and b. LDV measure- 
ments of the secondary flow velocities at one streamwise station 
(x ■ 229 mm) are presented and compared to the cross-wire measure- 
ments in Figs. 12a and b. LDV and cross-wire measurements of 
the turbulence quantities at one spanwise position (z ■ 13 mm) 
are compared in Fig. 13. 

The secondary flow velocities are qualitatively similar, 
although the LDV measurements indicate a somewhat higher w. 
The normal intensity w' 2 also seems slightly high. The higher 
w measurements are more likely caused by a slight misalignment 
of the beams relative to the tunnel axis rather than by remnants 
of the apparent coincidence problem (discussed above in Section 5) , 
since the latter would not affect the w measurements, v' 2 seems 
to agree very well whereas u'^ is a bit low, and since u was 
a bit high, this was probably a result of fringe biasing caused 
by too high a mixing frequency (as discussed above in Section 4.2) . 
(The mixing frequencies used for these measurements were 37.5 Mhz 
for the u channel, and 38 MHz for the v and w v channels.) However, 
the normal stress measurements agree to within 10%, and the 
shear stresses are consequently affected. u'v' is somewhat 
low (about 20%) whereas u'w 1 is low by almost a factor of two. 
The measurement of u'w 1 with the present system seems to be 
very sensitive. Thi3 is due to the fact the w v * u is of the 
same order as u * v cos 45° , so any small error in the measurement 
of these velocities or the off-axis angle can result in large 
errors in u'w'. v'w' in vortex affected flows is generally 
of the same order as u'w' and this seems to be the case with 
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the present measurements. With a cross-wire, v'w' has to be 
evaluated from measurements made in four different planes about 
the probe axis, and hence was not measured here. The comparisons 
clearly demonstrate the potential of the system in measuring 
detailed mean flow and turbulence quantities in three-dimensional 
flows. Work is in progress on optimizing the problems discu sed 
above so that the measurement accuracies may be improved. 

6. CONCLUDING REMARKS 

A 3-component LDV system, capable of measuring all three 
components of velocity simultaneously has been developed for 
use in low-speed three-dimensional flows. All the six components 
of Reynolds shear stress and higher order products of interest 
can hence be evaluated. The approach followed was to convert 
an existing 2-component system by splitting one of the colors 
to produce the third beam pair. The additional optical hardware 
required for this process is relatively minor. 

For the first time, three-component measurements made with 
an LDV system have been compared directly with those obtained 
with the cross-w^re technique, in a three-dimensional flow field 
where both techniques are expected to perform satisfactorily. 
The preliminary measurements are encouraging and work is in 
progress on improving the system accuracy. 
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APPENDIX 


SOFTWARE FOR THE HP 9845B DESK-TCP COMPUTER 


Complete listings of two programs written in 
ASSEMBLY languages are included in tnis appendix: 
aata acquisition, some on-line data reduction, and 
data on floppy disks; and "STAT" for complete off 
reduction from files written to disk. 


BASIC and 
M LDV" for 
storage of 
-line data 
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10 

20 

30 

40 

30 

60 

70 

80 

90 

100 

110 

120 

130 

140 

130 

160 

170 

180 

190 

200 

210 

220 

230 

240 

230 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

368 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

300 

510 

320 

330 

540 

350 

366 

370 

580 

390 

600 


REM PROGRAM LDV 

t PROGRAM TO ACQUIRE DATA FROM THREE-COMPONENT LDV SYSTEM 

• ~he program a.ik i For initialization data and calculates the 
1 calibration constants From them. It reads 3 channels oF raw 
' LV data From -.he LDV A^D Cl and unites them to a disk File 

1 togetner with the calibration constants. The program can reduce 

• the data and display real-time histograms oF the raw LV data 
' if desired. 

OPTION BASE 1 
ICOM 16600 
IDELETE ALL 
I ASSEMBLE Find_vel 
IASSEMBLE Data -.rans 
1 ASSEMBLE Draw 'hist 

COM INTEGER Dat a< 3 , 2608 > , D1 <2000, 6 > , Ns, Nn 'D1 is data buFfer 
COM REAL DFl,DF2,DF3 

COM REAL Theta, Nub , Num i xi , Num i x2, Numi x3 
COM REAL Su, Su, 3u, Suu, Svv, Sww, Suu, Suw, Svw 
DIM Date$t303,Fi 1 e#C6J , Name*! 4 3 , T 1 1 1 SC 80 3 
REAL Xpos, Ypos, Zpos, Pos<30, 3> 

INTEGER A, B, Rangel , Range2 
INTEGER Run,Dn,Nss,N 
REAL Phl,D,F,Ph2 

REAL Lam 1 , Lam2, Lam3, Db, Prwi dl , Pr 1 enl 
REAL Re, Ue 

REAL Prw l d2, Prw i d3 , Pr 1 en2, Prl en3 

REAL NFrl,NFr2,Nfr3 

REAL U,V,W 

Ns-2000 

DEG 

PRINT 

PRINT " ♦* << PROGRAM LDV : 3-COMPONENT VELOCITY DATA >> ** " 

PRINT 

PRINT "PROGRAM STRUCTURE” 

PRINT “ 1. INITIALIZE VARIABLES AND CALCULA TE PARAMETERS ". 

PRINT " 2. ACQUIRE DATA FROM A.-D" 

PRINT ” 3. WRITE TO FLOPPY DISC" 

PRINT 

l 

' ♦* CHECK HISTOGRAMS ** 

t 

Ans*»"N" 

INPUT "DO YOU WISH TO LOOK AT HISTOGRAMS ’ <YxN, DEFAULT N>",Ans* 
IF Ans*-“Y" THEN GOSUB Hist 

i 

' #* INITIALIZE RUN ** 

i 

PRINT ” ** INITIALIZATION ** " 

Ruri» 1 

PRINT " ENTER RUN PARAMETERS: M 
PRINT 

INPUT "Enter date and t i me: " , Dat e* 

INPUT "Enter 1-Line Name For Profile : " , T i t 1 S 

INPUT "No. of dtta samples per point <2000 maxm.) :",Ns 

i 

' BEAM SPACING J IS FIXED 

i 

D-.013 

' FOCAL LENGTH F IS FIXED 
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610 I 
620 F-.381 

630 I 

640 ' WAVELENGTHS OF 3 BEAMS ARE FIXED 

630 ' 

660 Lam1“4.88E-7 

670 Lam2*3. 143E-7 

680 La»3-3. 143E-7 

690 I 

700 ! REFERENCE VELOCITY SET TO ZERO 

710 ' 

720 Ue*0 
730 I 

740 ' CALCULATE HALF-ANGLES FROM BEAM SPACINGS AND FOCAL LENGTH 

730 ' 

760 Phl"ATN(D/2/F> 

7 7 0 Ph2»Phl*2 
780 ! 

799 ! CALCULATE FRINGE SPACINGS FROM WAVELENGTHS AND HALF-ANGLi 

800 • 

810 Df l-La»t/2/SIN(Phl ) 

820 Df2»L*»2/2/SIN(Phl > 

830 Df 3«LAB»3< / 2<'SIN<Phl ) 

840 I 

830 • BEAM DIAMETER Db IS FIXED 

860 • 

870 D6-1.20E-3 

880 • 

890 1 CALCULATE PP03E VOLUME WIDTH AND LENGTH 

900 ’ 

910 Prwidl-4*L*«l *F'PI/Db/COS<Phl > 

920 Pr 1 en 1 «4 *L am 1*F 'PI 'Db/SIN(Phl > 

930 Prw 1 d2*4*Lam2*F 'PI/Db/COS(Phl > 

940 Prlen2-4*l \m2*F 'PI/Db/SIN^hl > 

930 Pruid3»Prwid2 
960 Pr 1 en3*Pr 1 en2 
970 Nfr l*Prwi dl/Df 1 
980 Nfr2-«Prui d2/Df2 
990 Nfr3-Nfr2 
1000 • 

1010 I GET MORE RUN PARAMETERS 
1020 • 

1030 INPUT "Enter Br tgg shift frequency (MHz) : " , Nub 
1040 Kil1_u-3 

1050 INPUT "Is this i tuo-channel or three-channel run <2/3, default 3) ,B ,Kil' 
_u 

1*060 IF Kill _u*3 THEN GOTO 1110 

1070 INPUT "Enter mi<ing frequency (MHz, 2 nos. ) “ , Num i x 1 , Num i x2 
1080 Numix3«0 
1090 Theta-0 
1100 GOTO 1130 

1110 INPUT “Enter mi<ing frequency (MHz, 3 nos . 1 ” , Nun. l < 1 , Num i x2 , Num i x3 

1120 INPUT "Enter 3rd beam angle (degrs) Theta 

1130 INPUT “Enter tunnel reference voltage < vo I t s > " , Vref 

1140 PRINTER IS 0 

1 150 ■ 

1160 • PRINT HEADER 

1170 i 

1 180 PRINT Till* 

1190 PRINT "TEST DATE AND TIME :",Date* 

1200 PRINT "BEAM SPACINGS 'm> ■“,D 
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1210 

1220 

1230 

1240 

1239 

1260 

1270 

1200 

1290 

1300 

1310 

1320 

1330 

1340 

13S0 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1430 

1460 

1470 

1480 

1490 

1500 

1310 

1320 

1330 

1340 

1330 

1360 

1370 

1380 

1390 

1680 

1610 

1620 

1630 

1640 

1630 

1660 

1670 

1680 

1690 

1700 

1710 


PRINT “TOTAL ANGLE BETWEEN BERMS <Degrs) *",PH2 
PRINT “FRINGE SPRCINGS <m> , Df 1 , Df 2, Df3 

PRINT "PROBE WIDTHS <m> Prwi dl , Prwi d2, Prwi d3 
PRINT "PROBE LENGTHS On) - " , Pr 1 en 1 , Pr I «n2, Pr 1 en3 
PRINT "NO. OF FRINGES - " , Nf r 1 , Nf r2 , Nf r3 
PRINT "BRAGG SHIFT FREQUENCY <MHz> *",Nub 
PRINT "MIXING FREQUENCY (MHz) ■ " , Nui*t xl , Numi x2, Hum x3 
PRINT "THIRD BEAM SET ANGLE (Degrs) *“,Th«,-.i 
PRINT "RUN NUMBER - “,Run 

PRINT "TUNNEL REFERENCE VOLTAGE (volts) «",Vref 

PRINTER IS 16 

An*S*"N" 

INPUT "DO YOU WISH TO MAKE ANY CHANGES'’ (Y/N, DEFAULT N)", 
IF AnsS*"Y" THEN GOTO 540 

1 

! DATA ACQUISITION 
t 

PRINT 

PRINT " ** DATA ACQUISITION ** " 

REDIM D*t*<3,Ns> 

AnsS-"N“ 

INPUT "DO YOU WISH TO CHANGE FILE NAME ’ < Y/N, DEFAULT N)" 

IF AnsS*"N" THEi: GOTO 1490 

INPUT “Enter 4-digit f i I enamel", NameS 

PRINTER IS 0 

PRINT "4-digit -'ilename for profile NameS 
INPUT "Enter D 1 Number : “ , Dn 
PRINT "Disk Number »",Dn 
PRINTER IS 16 

INPUT "Enter X, Y, AND 2 1 oc at l ons : “ , Xpos , Ypos , Zpos 

l 

! TAKE DATA 

i 

PRINTER IS 0 

PRINT 

PRINT 

PRINT “POINT NUMBER IN PROFILE : ",Run 
PRINT “ X , Y , Z » Xpos, Ypos, Zpos 
PRINTER IS 16 
GOSUB Atod 

i 

i CALCULATE SAMPLE VELOCITIES 

i 

Ans*»"N" 

INPUT "DO YOU WISH TO OBTAIN ESTIMATES OF U AND V ’ (Y'N, 

IF AnsS»"N" THEN GOTO 1720 

RAD 

ICALL F i nd_vel 
GOTO 2420 

t 

I WRITE CALIBRATION CONSTANTS AND DATA TO DISK 


AnsS 


, AnsS 


DEFAULT N > “ 


1720 i 


1730 

1740 

1738 

1760 

1770 

1780 

1790 

1300 


GOSUB Df l 1 e 
AnsS-"Y" 

INPUT "DO YOU WISH 10 
IF AnsS*" N “ THEN GO CO 
Run*Run* 1 
AnsS* "N* 

INPUT "DO YOU WISH TO 
IF AnsS* “ Y " THEN GOTO 


TAKE ANOTHER POINT ’ <. Y/N , DEFAULT Y)",AnsS 
1820 


CHANGE ANY PARAMETERS ’ (Y-N, DEFAULT NV’.AnsS 
540 


, AnsS 
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1810 GOTO 1410 
1820 END 

1830 . ' M I M II I I I t II I I I I I I I I II I I I I I I I I I li t I I I I I t M I I I II I I I I I I II I 

1840 i *********** END OF MAIN PROGRAM LDV ********************* 

18S0 i 1111111)1111111111111111111111111111111111111111111111111 
I860 (Hod: I Subroutine for input from the LDV-A/D Cl 

1870 I Enter routine with Ns ■ no. simples 

1880 DISP “Press CONT to initiate data acquisition" 

1890 PAUSE 

1900 DISP -Acquiring Data" 

1910 RESET 10 

1920 CONTROL MASK 10; 1 

1930 WRITE 10 10, 3; 0 

1940 WRITE 10 10,3; 1 1 start handshake by setting CTL0 

1930 Nt *6*Ns 

i960 FOR 1*1 TO 3 

1970 Dumny*READBIN< 10) 

1980 NEXT I 

1990 REDIM D1 (Ns, 6 > 

2000 ' 

•2010 ENTER 10 WFHS Nt NOFORMAT; Dl(*> 'fast data acquisition 

2020 WRITE 10 10,3;0 

2030 PRINT 

2040 DIS^»Data acquisition complete" 

2030 PRINT 

2060 ICALL Data trans 

2070 RETURN 

2088 1 i i i i i i i ) m i i i i i it ii i i i i i i i i i i i i i t i i t ii i i i i i i i i i i i i i i i m i i i m i 

2890 Dfile: ' unite data file to floppy disc 

2100 PRINTER IS 16 
2110 PRINT 

2128 PRINT “** DATA FILE WRITE TO FLOPPY DISK **" 

2130 PRINT 

2140 PRINT “At this point be sure there is a floppy in drive 0 of" 

2150 PRINT "the 9895A with space for a file of 101, 256-byte records." 

2160 PRINT “Sje 

2170 PRINT 
2180 Ans*»“Y" 

2190 INPUT “DO YOU WISH TO WRITE THESE DATA TO DISK ? (Y/N, DEFAULT Y>",Ans* 
2208 IF Ans**"N“ THEN GOTO 2380 
2210 Fi 1 e**NameSLVAL* ( Run > 

2220 DISP "File “;File*;“ being written to disk" 

2230 MASS STORAGE IS ":H8,0,0" ' set floppy drive (9893A drive 8> as default 

2240 CREATE File*, 101 ' open file with 101 records 236 bytes each 

2250 ASSIGN File* TO #1 
2260 PRINT *1 ; Date* 

2270 PRINT *1 ; Tl 1 1 * 

2280 PRINT #1 ; Name* 

2290 PRINT ttl ; Dn 

2300 PRINT #1 ; Nub, N uim x 1 , Numi x2, Numi x3, Theta, Run 
2310 PRINT »lJVref,Ue,Df l,Df2,Df3 
2320 PRINT *1 ; Xpos, Ypos, ZpOS, Ns 
2330 MAT PRINT #i;Data 

2340 PRINT "*«*** File write completed ♦****" 

2330 ASSIGN * TO *1 ' close data file 

2360 MASS STORAGE IS ":H8,0,1" 1 reset program disk as mass storage 

2370 GOTO 2390 
2380 Run*Run-l 
2390 RETURN 

2400 ' l i i l i i i l i i l l i l l i i i i i i i i i i l i i l l l l l i i l i i i i i l l i i i i l i i I l i II l i l i l i 
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2410 Uguess: ' Estimatt U-co«ipon»nt of velocity 

2420 DISP "CALCULATING ESTIMATES OF U AND V VELOCITIES" 

2430 > 

2440 • 

2430 IF KMl_u-3 THEN GOTO 2490 

2460 ORTA 0,0,0 

2470 RESTORE 2460 

2480 READ Su,Swu,Suw 

2490 Ubar-SuxNs 

2300 Vbar-SvxNs 

2310 Wb*r«Sw''N» 

2320 Upr l 2»Suu''Ns-Ub »r*Ubar 
2530 Vpr i 2*SvvxNs-Vb »r*Vbar 
2340 Wprf 2»SwwxNs-Wb »r«Wbar 
2330 Uvbar»SuvxNs-Ub areVbar 
2560 Uwbar»SuwxNs-Ub »r*Wb*r 
2570 Vwbar“Svw/'Ns-Vb tr*Wbar 
2330 PRINTER IS 0 

2590 PRINT "ESTIMATE OF QUANTITIES FROM SAMPLES : M , N» 

2600 PRINT "Ubar -"jUbar 
2610 PRINT "Vbar -".‘/bar 
2620 PRINT "Ubar *",Ubar 
2630 PRINT "Uprt 2 -*,Upri2 
2640 PRINT " Vpr i 2 -",Vpri2 
2650 PRINT "Wpr i 2 »",Wpri2 
2660 PRINT "Uvbar Uvbar 
2670 PRINT "Uubar Uubar 
2630 PRINT “ Vwbar »",Vubar 
2690 GOTO 1700 

2700 ISOURCE NAM Find vel 

2710 ISOURCE i 

2720 ISOURCE 1 This subroutine converts raw data counts into l nst ant aneous 
2730 ISOURCE 1 velocities, then sums several different products of the 
2740 ISOURCE 1 velocity components. All input and output data is passed 
2750 ISOURCE 1 through the COMMON storage area. The inputs are the raw 
2760 ISOURCE 1 data trray (Arrayd), the Bragg shift frequency (Nub), the 
2770 ISOURCE 1 mixing frequencies (Numixl, Numix2, Num x3), the fringe 
2730 ISOURCE 1 spacings (Dfl, Df2, Df3>, the crossing angle of the third 
2790 ISOURCE 1 beam ITheta), and the number of samples in a data point (Ns). 
2800 ISOURCE 1 The outputs are the summations of various products of the 
2810 ISOURCE 1 velocity components, including U, V, U, U*U, V#V, W*W, U*V, 
2820 ISOURCE ■ U*W, md V*U. 

2330 ISOURCE • 

2840 ISOURCE EXT Get_value 1 Declare subroutines stored 

2850 ISOURCE EXT Get~info 1 outside of the main program. 

2860 ISOURCE EXT Get_element 

2870 ISOURCE EXT Int_to_rel 

2880 ISOURCE EXT Rel_math 

2890 ISOURCE EXT Put value 

2900 ISOURCE • 

2910 ISOURCE COM 

2920 ISOURCE Data par: INT <*) 1 Declare common variables. 

2930 ISOURCE D1 par: INT <*) 

2940 ISOURCE Ns~par: INT 

2930 ISOURCE Nn~par: INT 

2960 ISOURCE Dfl_par! REL 

2970 ISOURCE Df2_par REL 

2980 ISOURCE Df 3 par: REL 

2990 ISOURCE Theta ptr: REL 

3000 ISOURCE Nub_par : REL 
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3910 

ISOURCE 

Nmi xl_par: 

REL 



3020 

ISOURCE 

Nm i x2_p »r : 

REL 



3030 

I SOURCE 

Nmtx3_psr: 

REL 



3040 

ISOURCE 

Su_par : 

REL 



3050 

ISOURCE 

So_par: 

REL 



3060 

ISOURCE 

Sw par: 

REL 



3070 

ISOURCE 

Suu_par : 

REL 



3080 

ISOURCE 

S uu_p*r : 

REL 



3090 

ISOURCE 

Suu_par : 

REL 



3100 

ISOURCE 

Suo_par : 

REL 



3110 

ISOURCE 

Suu_par : 

REL 



3120 

ISOURCE 

Sou pv: 

REL 



3130 

ISOURCE 

! 




3140 

ISOURCE 

Arr ayd: 

BSS 

39 

1 Reseroe space for dAtA Array 

3150 

ISOURCi. 

El ementd: 

EQU 

Arrayd+16 ' descriptor. 

3160 

ISOURCE. 

Array l ! 

BSS 

4096 

1 Reseroe space for lookup tables used 

3170 

ISOURCE 

P**r av2: 

BSS 

4096 

! for count to velocity conversion. 

3180 

ISOURCE 

Array3: 

BSS 

4096 


3190 

ISOURCE 

Ns: 

BSS 

1 

• Reseroe space for oanous input 

3200 

ISOURCE 

Dfi: 

BSS 

4 

1 and output variables. 

3210 

ISOURCE 

Df2: 

BSS 

4 


3220 

ISOURCE 

Df3: 

BSS 

4 


3230 

ISOURCE 

Th#t a: 

BSS 

4 


3240 

ISOURCE 

Nub: 

BSS 

4 


3250 

ISOURCE 

Nu«i xl : 

BSS 

4 


3260 

I SOURCE 

Nunn x2: 

BSS 

4 


3270 

ISOURCE 

Nunn x3: 

BSS 

4 


3280 

ISOURCE 

Su: 

BSS 

4 


3290 

ISOURCE 

So: 

BSS 

4 


3300 

ISOURCE 

Su: 

BSS 

4 


3310 

ISOURCE 

Suu: 

BSS 

4 


3320 

ISOURCE 

Sou: 

BSS 

4 


3330 

ISOURCE 

Suu: 

BSS 

4 


3340 

ISOURCE 

Suo: 

BSS 

4 


3350 

ISOURCE 

Suu: 

BSS 

4 


3360 

ISOURCE 

Sou: 

BSS 

4 


3370 

ISOURCE 

Count : 

BSS 

1 ' 

Count and I are general purpose index vari- 

3388 

ISOURCE 

I: 

BSS 

1 ' 

ables. Count is usually 0, 1, or 2, to denote 

3390 

ISOURCE 

Check : 

BS~ 

1 ' 

whether U,V, or U is being calculated. 

3400 

ISOURCE 

Int : 

BSS 

1 ' 

Int, Address, and Offset are all general 

3410 

ISOURCE 

Address : 

BSS 

1 ' 

purpose storage areas. 

3420 

ISOURCE 

Offset : 

BSS 

1 


3430 

ISOURCE 

Rl: 

BSS 

4 ' 

Rl, R2, and R3 are the count-to- 

3440 

ISOURCE 

R2: 

BSS 

4 ' 

frequeney conversion factors. 

3450 

ISOURCE 

R3: 

BSS 

4 


3460 

ISOURCE 

Xvar : 

BSS 

4 ' 

Xoar and Yoar are general purpose real 

3470 

ISOURCE 

YoAr*. 

BSS 

4 ' 

number storage areas. 

3480 

ISOURCE 

u: 

BSS 

4 ' 

reseroe space for instantaneous velocity 

3490 

ISOURCE 

v: 

BSS 

4 ' 

components. 

350£< 

ISOURCE 

w: 

BSS 

4 


3510 

ISOURCE 

uu: 

BSS 

4 


3520 

ISOURCE 

Uo: 

BSS 

4 


3530 

ISOURCE 

Cos: 

BSS 

4 ' 

Cos and Sin are the cosine and sine of Thet* 

3540 

ISOURCE 

Sin: 

BSS 

4 


3558 

ISOURCE 

Rad: 

DAT 

5. 729578E 1 

3560 

ISOURCE 

Mill: 

DAI 

1.E6 


3570 

ISOURCE 

One: 

DAT 

1 . 


3583 

ISOURCE 

Zero: 

DAT 

0. 


3590 

ISOURCE 

i 




3600 

ISOURCE 


SUB 
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1 Get number of samples 


3610 

ISOURCE 

Find 

ve 1 : 

LOR 

-Ns 

3620 

ISOURCE 




LOB 

-Ns_par 

3630 

ISOURCE 




JSM 

Get_val ue 

3640 

ISOURCE 




LOR 

-flrrayd 

3650 

ISOURCE 




LOB 

-Oat a_par 

3660 

ISOURCE 




JSM 

Get info 

3670 

ISOUPCE 




LOR 

•Of l 

3680 

ISOURCE 




LOB 

-Of l_par 

3690 

ISOURCE 




JSM 

Get value 

3700 

ISOURCE 




LOR 

-Of 2 

3710 

ISOURCE 




LOB 

-Of 2_par 

3720 

ISOURCE 




JSM 

Get value 

3730 

ISOURCE 




LOR 

-Of 3 

3740 

ISOURCE 




LOB 

•Df3_par 

3730 

ISOURCE 




JSM 

Get_val ue 

3760 

ISOURCE 




LOR 

-Thet a 

3770 

ISOURCE 




LOB 

■Theta_par 

3780 

ISOURCE 




JSM 

Get_val ue 

3790 

ISOURCE 




LOR 

-Nub 

3800 

ISOURCE 




LOB 

-Nub_par 

3810 

ISOURCE 




JSM 

Get_val ue 

3820 

ISOURCE 




LOR 

■Numi xl 

3830 

ISOURCE 




LOB 

-Nmi xl_par 

3840 

ISOURCE 




JSM 

Get_val ue 

3830 

ISOURCE 




LOR 

-Num i x2 

3860 

ISOURCE 




LOB 

-Nmi x2_par 

3870 

ISOURCE 




JSM 

Get_val ue 

3880 

ISOURCE 




LOR 

-Num i x3 

3890 

ISOURCE 




LOB 

-Nmi x3_par 

3900 

ISOURCE 




JSM 

Get_val ue 

3910 

ISOURCE 

1 





3920 

ISOURCE 

1 

The 

1 oop 

headed by Get f 

3930 

ISOURCE 

1 

the 

count 

-t o- 

■frequency c 

3940 

ISOURCE 

t 

the 

range) for U, V, and 

3950 

ISOURCE 

1 

the 

v ir i ab 1 e 

" Count ", th 

3960 

ISOURCE 

1 

the 

s im* 

for 

U, V, and W 

3970 

ISOURCE 

1 





3980 

I SOURCE 




LOR 

-0 

3990 

ISOURCE 




STfl 

Count 

4000 

I SOURCE 

Get freg: 

LOR 

Cou 

4010 

ISOURCE 




LOB 

Ns 

4020 

ISOURCE 




MPY 


4030 

ISOURCE 




STR 

El ementd 

4040 

ISOURCE 




LOR 

-Int 

4050 

ISOURCE 




LOB 

-flrrayd 

4060 

ISOURCE 




JSM 

Get_e 1 ement 

4070 

ISOURCE 




LDfl 

Int 

4080 

ISOURCE 




LOB 

-13360 

4090 

ISOURCE 




RND 

B 

4100 

ISOURCE 




SRR 

10 

4110 

ISOURCE 




TCR 


4120 

ISOURCE 




LOB 

-15 

4130 

ISOURCE 




ADR 

B 

4140 

ISOURCE 




LOB 

-1 

4130 

ISOURCE 




SZR 

Loopend 

4160 

ISOURCE 

Loop: 


SBL 

1 

4170 

ISOURCE 




DSZ 

R 

4180 

ISOURCE 




JMP 

Loop 

4190 

ISOURCE 

Loopend : 

STB 

I nt 

4200 

ISOURCE 




LDfl 

-Int 


Get parameter* of data array. 


! Get input parameters. 


Whenever a loop is controlled by 


• Get the first word of the column of 

• the data array which contains the 

1 velocity component for which we want 
' to get the range. 


• Mask and rotate to get the four 

• bits containing the range. 


' Subtract from 13 to get the 
1 actual range. 


1 Use the range to find the power 
i of two needed for the divisor. 


1 Convert the power of two into a 


* 



/ 



4218 

ISOURCE 

STfl 

Gper_l 

1 real number. 

4228 

ISOURCE 

LOR 

■ Yvar 


4238 

ISOURCE 

STfl 

Result 


4240 

ISOURCE 

JSM 

Int to ret 


42S0 

ISOURCE 

LDR 

--3?2E4 

1 Divide 3.2E4 by the appropriate 

4260 

ISOURCE 

LDB 

-Xvar 

1 power of two, using BCD math. 

4270 

ISOURCE 

XFR 

4 


4280 

ISOURCE 

STB 

Oper_l 


4290 

ISOURCE 

LOR 

■Yu*!* 


4300 

ISOURCE 

STfl 

0per_2 


4310 

ISOURCE 

LOR 

Count 

1 Decide whether to put the result 

4320 

ISOURCE 

LOB 

■4 

1 in Rl, R2, or R3, depending on 

4330 

ISOURCE 

MPY 


' Count. 

4340 

ISOURCE 

ADR 

*R1 


43S0 

ISOURCE 

STfl 

Result 


4360 

ISOURCE 

LOR 

■2 


4370 

ISOURCE 

LOB 

-147155B 

• Now, finally, call the utility to 

4380 

ISOURCE 

JSM 

Re1_math 

* perform the division. 

4390 

ISOURCE 

ISZ 

Count 


4400 

ISOURCE 

LOR 

■ 3 

' Increment and check Count so as 

4410 

ISOURCE 

CPA 

Count 

1 to follow the loop three times. 

4420 

ISOURCE 

JMP 

*+2 


4430 

ISOURCE 

JMP 

Get_f neq 


4440 

ISOURCE ' 




4450 

ISOURCE 

LOR 

-flrrayl 

• Zero out the entire count-to- 

4460 

ISOURCE 

LOB 

-768 

1 velocity conversion table so that 

4470 

ISOURCE Continue: 

CLR 

16 

' it must be recalculated for each 

4480 

ISOURCE 

ADR 

-16 

1 point. (This must be done if the 

4490 

ISOURCE 

DSZ 

B 

1 mixing frequencies or ranges are 

4500 

ISOURCE 

JMP 

Cont i nue 

1 changed between counts. > 

4510 

ISOURCE ' 




4520 

ISOURCE 

LOR 

-Su 

' Set initial values of Su, Sv, 

4530 

ISOURCE 

LDB 

-9 

1 Su, Suu, etc. to zero. 

4540 

ISOURCE C 1 ear : 

CLR 

4 


4550 

ISOURCE 

ADR 

-4 


4560 

ISOURCE 

OSZ 

B 


4570 

ISOURCE 

JMP 

Clear 


4580 

ISOURCE ' 




4590 

ISOURCE 

LOR 

-Theta 

1 Convert Theta from degrees to 

4600 

ISOURCE 

STfl 

Oper_l 

• radians using the Rel_math 

4610 

ISOURCE 

LOR 

-Rad 

■ utility. 

4620 

ISOURCE 

STfl 

0per_2 


4630 

ISOURCE 

LOR 

-Xuar 


s640 

ISOURCE 

STfl 

Result 


4650 

ISOURCE 

LOR 

-2 


4660 

ISOURCE 

LOB 

-147155B 


4670 

ISOURCE 

JSM 

Rel_math 


4688 

ISOURCE ' 




4690 

ISOURCE 

LOR 

-Xuar 

• Find the sine and cosine of 

4700 

ISOURCE 

STfl 

Oper_t 

1 Theta, and store them in the 

4710 

ISOURCE 

LOR 

-Si n 

1 locations Sin and Cos, respect i vel y. 

4720 

ISOURCE 

STR 

Result 


4730 

ISOURCE 

LOR 

-1 

1 

4740 

ISOURCE 

LDB 

- 34213 B 


4750 

ISOURCE 

JSM 

Rel_math 


4760 

ISOURCE 

LOR 

-Cos 

* 

4770 

ISOURCE 

STfl 

Resu 1 t 


4780 

ISOURCE 

LOR 

-1 


4790 

ISOURCE 

LOB 

-34224B 


4800 

ISOURCE 

JSM 

Rel math 
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ISOURCE 
I SOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 


Thv I' op defined by Get_int calculate* intermediate values 
used i i converting counts to velocities. These values are the 
same -cm all samples in a point, so they can be calculated 
separately. The loop calculates <Nub _ NumixN> and <Mi 1 1 *DfN> , 
where N is 1, 2, and 3. 


Get int 


LDR *3 
STR Count 
LDR Count 
RDR —1 
SRL 2 
STR Offset 
LDR ■Nub 
STR Oper_l 
LDR ■Numixl 
RDR Offset 
STR 0per_2 
STR Result 
LDR -2 

LDB “146717B 
JSM Rel_math 
LDR ■Df 1 
RDR Offset 
STR Oper_l 
STR Result 
LDR “Mi 1 1 
STR Open 2 
LDR -2 

LDB “147637B 
JSM Re1_math 
DSZ Count 
JMP Get int 


• Find Numi xN“Nub-Numi xN. 


Find DfN-Mi 1 1 *DfN. 


The loop Begin 
and W. Each ti 
raw d it urn) is 
The velocity i 
this is the fi 
first timt any 
corre tpondi ng 
found in Get_i 
count is found 
oc l ty is 1 ooke 


i s perf ormed 
me through th 
taken from it 
s stored in U 
rst, second, 
particular c 
to it is cal c 
nt, and the v 
again in the 
d up rather t 


three times, once for U, V, 
e data array is read, a count (a 
. and is converted into a velocity. 

, V, or U depending on whether 
or third iteration of the loop. The 
ount is encountered, the velocity 
ulated using the intermediate values 
elocity is stored in a table. If that 
data array, the correspondi ng vel- 
han being calculated again. 


Begin! 


LDR “0 
STR Count 
LDR “5 
STR Check 
LDR Ns 
•ITR I 

l DR Count 
STR B 
3BL 2 
3TB Offset 
LDB Ns 
MPY 
RDR I 
RDR —1 


• Figure out which element of the 

• data array we want to pick up. 


I 
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l 


T 

i 

f 


r 
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'y. 


i 
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L 
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3410 

I SOURCE 

STA 

El ementd 

3420 

ISOURCE 

LDA 

■ Int 

5430 

ISOURCE 

LDB 

■Arrayd 

5440 

ISOURCE 

JSM 

Get_e 1 ement 

3450 

ISOURCE 

LDA 

Int 

3460 

ISOURCE 

LDB 

>1023 

3470 

ISOURCE 

AND 

B 

5480 

ISOURCE 

LDB 

Count 

3490 

ISOURCE 

TCB 


3500 

ISOURCE 

ADB 

Check 

3510 

ISOURCE 

SZB 

Strai ght 

5520 

ISOURCE 

TCA 


5330 

ISOURCE 

ADA 

-1024 

5540 

ISOURCE Straight: 

STA 

Int 

3330 

ISOURCE 

LDA 

Count 

5560 

ISOURCE 

LDB 

■4096 

5570 

ISOURCE 

MPY 


5580 

ISOURCE 

LDB 

Int 

5590 

ISOURCE 

ADB 

— 1 

3600 

ISOURCE 

SBL 

2 

5610 

ISOURCE 

ADA 

B 

3620 

ISOURCE 

ADA 

•l 

3630 

ISOURCE 

ADA 

■Array 1 

3640 

ISOURCE 

STA 

Address 

5650 

ISOURCE 

LDA 

Address, I 

3660 

ISOURCE 

SZA 

Cal cu 1 at e 

3670 

ISOURCE 

JMP 

Over 

5680 

ISOURCE CAlculAte: 

LDA 

■ Int 

3690 

ISOURCE 

STA 

Oper_l 

5700 

ISOURCE 

LDA 

■ Yvar 

3710 

ISOURCE 

STA 

Result 

3720 

ISOURCE 

JSM 

Int_t o_rel 

3730 

ISOURCE 

STB 

Oper 2~ 

5740 

ISOURCE 

LDA 

-R 1 

5750 

ISOURCE 

ADA 

Offset 

3760 

ISOURCE 

STA 

Oper_l 

5770 

ISOURCE 

LDA 

■Xvar 

3780 

ISOURCE 

STA 

Result 

5790 

ISOURCE 

LDA 

-2 

5800 

isourcE 

LDB 

■ '.47135B 

5810 

ISOURCE 

JSM 

Rt l_math 

5820 

ISOURCE i 



3830 

ISOURCE 

LDA 

■Numi xl 

3840 

ISOURCE 

ADA 

Offset 

3830 

ISOURCE 

STA 

Oper_l 

5860 

ISOURCE 

LDA 

■Xvar 

5870 

ISOURCE 

STA 

0per_2 

5880 

ISOURCE 

LDA 

■ Yvar 

5890 

ISOURCE 

STA 

Result 

3900 

ISOURCE 

LDA 

■2 

3910 

ISOURCE 

LDB 

■146717B 

5920 

ISOURCE 

JSM 

Rel_ math 

3930 

ISOURCE • 



3940 

ISOURCE 

LDB 

Count 

5950 

ISOURCE 

ADB 

— 1 

3960 

ISOURCE 

SBM 

*+2 

3970 

ISOURCE 

JMP 

Sames i gn 

3980 

ISOURCE 

LDA 

■Zero 

3990 

ISOURCE 

STA 

Oper_l 

6000 

ISOURCE 

LDA 

ResuTt 


! Get a r*u datum from t he d At a array. 


1 Strip off the ftrst six 
1 bits of the raw d At a word. 
' See if Count » Check. 


I If true, use the stripped dAtA 
! word as An index. If not, use 
I 1024 minus the dAtA word. 

1 Store the count we hAve gotten. 

■ Now, use Count to find out 
! which lookup tAble array 
' we wAnt to use. And use the 
I count we got from the dAtA ArrAy 
I to find exActly where in the tAble 
! we WAnt to go. 


• If thAt tAble entry is zero, 
! cAlculAte a velocity for it. 


' Convert the count to a reAl number. 


' Divide the rang* we found 
1 ear 1 i«r by the count to get a 
1 frequency. 


1 Find <Nub-Numi xN>-FrequencyN. 


• If we Are cAlculAting U, reverse 
1 the sign of <<Nub-NumixN)-Frequ*ncyN'' 
1 so as to reverse the sign of U. 

1 Leave V and W Alone. 
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STR 0per_2 
LDR *Xvar 
STR Result 
LDR -2 

LOB *1467176 
JSM Re 1 math 
NOP 


LDR Result 
STR Oper 1 
LDR -Dfl“ 

RDR Offset 
STR 0per_2 
LDR Rddress 
RDR —1 
STR Result 
LDR *2 

LDB * 1 47037B 
JSM Rel math 


Find Vel oc 1 ty*< <Nub-Numl xN> 
-FrequencyN>*<Mi 1 1 *Df N> 
end store in * piece in the 
lookup teble correspondi ng tc 
the dete count. 


LDR Rddress 
RDR *-l 
LDB *U 
RDB Offset 
XFR 4 


1 Now trensfer the velocity from 
I the lookup teble to U, V, or 14, 
* es eppropnete. 


ISZ Count 
LDR ®-3 
RDR Count 
NOP 

SZR *+2 
JMP Begin 
STR Count 


1 Have U,V, end 14 ell 
1 been calculated 1 * 

1 If not, go beck again. 

1 If so, set Count * 8. 


Now we convert the w we have obtained (which is measured at 
an angle Theta to the V-axis> to the 14 we want (which should 
be meisured at an angle of 98 degre_sa_ijO — the V-raxi sH 
Thus -'ind 14* <l4v-V*Cos<Theta) )/Si n<Theta) . 


*V 

Oper_l 

■ Cos 
0per_2 
■Xver 
Resul t 
-2 

■147037B 
Re 1 _mat h 
Resul t 
Oper 2 

-W 

Oper_l 

■ Yvar 
Result 
-2 

-146717B 
Re 1 _mat h 
Resul t 
Oper_l 
*S l n 


/ 


i 



/ 


* 


6610 
6620 
6630 
6640 
6650 
6660 
6670 
6680 
6690 
6700 
6710 
6720 
6730 
6740 
6750 
6760 
6770 
6780 
6790 
6800 
6810 
6820 
6830 
6840 
6850 
6860 
6873 
6880 
6890 
6900 
6910 
6926 
6930 
6940 
6950 
6960 
6970 
6980 
6990 
7000 
7010 
7020 
7030 
7040 
7050 
7060 
707 a 
7080 
7090 
7100 
7110 
7120 
7130 
7140 
7150 
7160 
7170 
7180 
7190 
7200 


ISOURCE 

I SO'JRCE 

I SOURCE 

I SOURCE 

ISOURCE 

ISOURCE 

ISOURCE I 

ISOURCE ! 

ISOURCE • 

ISOURCE ' 

ISOURCE ' 

ISOURCE • 

ISOURCE I 

ISOURCE I 

ISOURCE • 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE ! 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOUPCE I 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE ' 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE ' 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE ! 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE ' 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE ' 

ISOURCE 

ISOURCE 

I SOURCE 

ISOURCE ' 


STA Oper 2 
LDA *14 
STA Result 
LDA *2 

LDB -147155B 
JSM Rel_m*th 

New tike running sums of U, V, 14, And several produces 
of these velocities. The sums Are taken using the utility 
"Add". The sums Are CAlculAted in An unusuAl sequence in 
order to reduce the number of steps needed to CAlculAte 
them. ProgrAm steps which produce sums not considered 
necessary have been included, but these steps Are preceded 
by the " 1 " comment »*rk«r. 

LDA *Su 1 Find Su*Su+U 

LDB -U 
JSM Add 

LDP Oper 2 I Find U*U 

STA Oper 1 

LDA -Uu 

STA Result 

LDA *2 

LDB -147037B 
JSM Rel_mAth 

LDA -Suu ' Find Su-Su+<U*U> 

LDB Result 
JSM Add 

LDA -Sv 1 Find Sv*Sv*'' 

LDB *V 
JSM Add 

LDA -U ' Find U#V 

STA Oper 1 

LDA *Uv 

STA Result 

LDA -2 

LDB *147037B 
JSM Rel_mAth 

LDA *Suv 1 Find Suv*Suv+<U*V> 

LDB Result 
JSM Add 

LDA *V ! Find V*V 

STA Qper_l 

STA 0per_2 

LDA *Xvar 

STA Result 

LDA -2 

LDB *147037B 
JSM Rel_math 

LDA *Svu ' Find Suv*Svv+< V*V> 

LDB Result 
JSM Add 





7210 

ISOURCE 

LDfl 

-Su 

' Find 

Sw"Sw+W 

7220 

ISOURCE 

LOB 

-M 



7230 

ISOURCE 

JSM 

Add 



7240 

ISOURCE ' 





7250 

ISOURCE 

LOR 

»U 

• Ft na 

u*w 

7260 

ISOURCE 

STfl 

Qper_i 



7270 

ISOURCE 

LOR 

»Xvar 



7280 

ISOURCE 

STfl 

Result 



7290 

ISOURCE 

LOR 

-2 



7300 

ISOURCE 

LOB 

■147037B 



r'310 

ISOURCE 

JSM 

Re1_math 



7320 

ISOURCE ' 





7330 

ISOURCE 

LOR 

•Suu 

• Find 

Suw*Suw+<U*W> 

7340 

ISOURCE 

LOB 

Result 



7350 

ISOURCE 

JSM 

Add 



7360 

ISOURCE ' 





7370 

ISOURCE 

LOR 

*U 

• Find 

V*W 

7380 

ISOURCE 

STfl 

Oper 1 



7390 

ISOURCE 

LOR 

■ V 



7400 

ISOURCE 

STfl 

Gper_2 



7410 

ISOURCE 

LOR 

■ Xuar 



7420 

ISOURCE 

STfl 

Result 



7430 

ISOURCE 

LOR 

-2 



7440 

ISOURCE 

LOB 

* 1 47037B 



7450 

ISOURCE 

JSM 

Rel_math 



7460 

ISOURCE ' 





7470 

ISOURCE 

LOR 

*Svw 

* Find 

Svw«Svu+< V*W> 

7480 

ISOURCE 

LOB 

Result 



7490 

ISOURCE 

JSM 

Add 



7500 

ISOURCE ' 





7510 

ISOURCE 

LOR 

-W 

' Find 

W*U 

7520 

ISOURCE 

STfl 

Oper_l 



7330 

ISOURCE 

STfl 

0per_2 



7540 

ISOURCE 

LOR 

■Xvar 



7550 

ISOURCE 

STfl 

Result 



7560 

ISOURCE 

LDfl 

■2 



7570 

ISOURCE 

l 

■ 147037B 



7580 

ISOURCE 

ISM 

R*!_jnath 



7590 

ISOURCE 1 





7600 

ISOURCE 

-or, 


• Find 

Sww*Swu+<W*W> 

7610 

ISOURCE 

LOB 

Resul t 



7620 

ISOURCE 

JSM 

Add 



7630 

ISOURCE ' 





7640 

ISOURCE 

DSZ 

I 

Cont 

nue to calculate running 

7650 

ISOURCE 

JMP 

Begin 

* sum* 

vn'.U out of samples. 

7660 

ISOURCE ' 






7670 ISQURCE 1 Place the finished sum in the COMMON region so that 
7680 ISOURCE 1 the BASIC program has access to them, and then return to 


7690 

ISOURCE 

the BASIC program. 

7^00 

ISOURCE 



7710 

ISOURCE 

LOR 

*Su 

77 '. 0 

ISOURCE 

LOB 

■Su_par 

7730 

ISOURCE 

JSM 

Put value 

7740 

ISOURCE 

LDfl 

-Sv 

7750 

ISOURCE 

LDB 

=Sv_par 

7760 

ISOURCE 

JSM 

Put_val ue 

7770 

ISOURCE 

LDfl 

■ Sw 

7780 

ISOURCE 

LDB 

*Sw_par 

7790 

ISOURCE 

JSM 

Put_val ue 

7800 

ISOURCE 

LDfl 

*Suu 


7810 

ISOURCE 


LOB 

»Suu_par 



7820 

1 SOURCE 


JSM 

Put_val ue 



7830 

ISOURCE 


LDR 

>Svv 



7840 

ISOURCE 


LOB 

*Svv_par 



7850 

ISOURCE 


JSM 

Put_val ue 



7860 

ISOURCE 


LOR 

*Suu 



7870 

ISOURCE 


LOB 

»Sww_par 



7880 

ISOURCE 


JSM 

Put_val ue 



7890 

ISOURCE 


LOR 

■Suv 



7900 

ISOURCE 


LOB 

=*Suv_par 



7910 

ISOURCE 


JSM 

Put_val ue 



7920 

ISOURCE 


LOR 

*Suw 



7930 

ISOURCE 


LOB 

*Suw_par 



7940 

ISOURCE 


JSM 

Put _val ue 



7950 

ISOURCE 


LOR 

=>Svw 



7960 

ISOURCE 


LOB 

=Svu_par 



7970 

ISOURCE 


JSM 

Put value 



7980 

ISOURCE 


RET 

1 



7990 

ISOURCE 

1 





8000 

ISOURCE 

• The i 

1 1 1 y " 

Add" is used 

to add up 

the running sums. 

8018 

ISOURCE 

i 





8820 

I SOURCE 

Add: 

ISZ 

Ut 1 count 



8030 

ISOURCE 


STR 

Oper_l 



8040 

ISOURCE 


STB 

0per_2 



8050 

ISOURCE 


STR 

Resul t 



8060 

ISOURCE 


LOR 

*2 



8070 

ISOURCE 


LOB 

■ 1 4672 1 B 



8080 

ISOURCE 


JSM 

Rel_math 



8090 

ISOURCE 


OSZ 

Ut 1 count 



8100 

ISOURCE 


RET 

1 



8110 

ISOURCE 


JSM 

Ut 1 end 



8120 

ISOURCE 

i 





8130 

ISOURCE 


LiT 

200 



8140 

ISOURCE 


END 

F i nd_ve 1 



8150 

ISOURCE 

i 





8*50 

ISOURCE 


NRM 

Oat a_t rjns 



8170 

ISOURCE 

i 





"ST30 - 

isour:e 

1 Subroutine to transfer raw 

data from input array 

8190 

ISOURCE 

1 to storage array. 



8200 

I SOURCE 

i 





8210 

ISOUFCE 


EXT 

Get_’nfo 1 

Dec 1 are 

subroutines stored 

8222 

ISOURCE 


EXT 

Get_va1ue 1 

out s i de 

of the main program 

8230 

ISOURCE 


tXT 

Tet _e 1 ement 



8240 

ISOURCE 


EXT 

Hut ^el ement 



8250 

ISOURCE 

i 





8260 

ISOURCE 


COM 




8270 

ISOURCE 

Data_par: 

INT 

(*) • 

Dec 1 are 

common variables. 

82£ 0 

ISOURCE 

Dl_par : 

INT 

<*> 



8290 

ISOURCE 

Ns_par : 

INT 




8300 

ISOURCE 

i 





8310 

ISOURCE 

Rrrayd: 

BSS 

39 



8320 

ISOURCE 

El ementd: 

EQU 

firrayd+16 



8330 

ISOURCE 

Array l : 

BSS 

39 



83-te 

ISOURCE 

E 1 ement l : 

EQU 

flrrayl+16 



8350 

ISOURCE 

Ns: 

BSS 

1 



8 ''60 

ISOURCE 

Count 1 : 

BSS 

1 



8-70 

ISOURCE 

Count 2: 

BSS 

1 



8380 

ISOURCE 

I : 

BSS 

1 



8390 

ISOURCE 

Int: 

BSS 

1 



8400 

ISOURCE 

1 
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*1 

*?i 


I? 


■j 


8410 

isource 


SUB 


8429 

ISOURCE 

Dat A_tr ins 

: LDA -Arrayd 

8430 

ISOURCE 


LOB 

■ Dat a_par 

8440 

ISOURCE 


JSM 

Get_i nfo 

8450 

ISOURCE 


LDA 

■Arr ay 1 

8460 

ISOURCE 


LOB 

*Dl_par 

8470 

ISOURCE 


JSM 

Get info 

8480 

ISOURCE 


LDA 

■Ns“ 

8490 

ISOURCE 


LDB 

■Ns_par 

8300 

ISOURCE 


JSM 

Get_ual ue 

8510 

ISOURCE 

i 



8320 

ISOURCE 


LDA 

-0 

8530 

ISOURCE 


STA 

I 

8340 

ISOURCE 

St art l : 

LDA 

-0 

8550 

ISOURCE 


STA 

Count 1 

8360 

ISOURCE 


LDA 

-2 

8370 

ISOURCE 


STA 

Count 2 

8580 

ISOURCE 

i 



8590 

ISOURCE 

St art2: 

LDA 

I 

8680 

ISOURCE 


LDB 

*6 

8610 

ISOURCE 


MPY 


8620 

ISOURCE 


ADA 

Count 2 

8630 

ISOURCE 


STA 

E 1 ement 1 

8640 

ISOURCE 


LDA 

-Int 

8650 

ISOURCE 


LDB 

■Arr ay 1 

8660 

ISOURCE 


JSM 

Get_e 1 ement 

8670 

ISOURCE 

i 



8680 

ISOURCE 


LDA 

Count 1 

8690 

ISOURCE 


LDB 

Ns 

8700 

ISOURCE 


MPY 


8710 

ISOURCE 


ADA 

I 

8720 

ISOURCE 


STA 

El ementd 

8730 

ISOURCE 


LDA 

■ Int 

8740 

ISOURCE 


LDB 

■Arr ayd 

8750 

ISOURCE 


JSM 

Put_e 1 ement 

8760 

ISOURCE 


ISZ 

Count 1 

8770 

ISOURCE 


ISZ 

Count 2 

8780 

ISOURCE 


LDA 

--3 

8790 

ISOURCE 


ADA 

Count l 

8800 

ISOURCE 


SZA 

*+2 

8810 

ISOURCE 


JMP 

St art 2 

8820 

ISOURCE 

i 



3830 

ISOURCE 


ISZ 

I 

8840 

ISOURCE 


LDA 

I 

8830 

ISOURCE 


TCA 


8860 

ISOURCE 


ADA 

Ns 

8870 

ISOURCE 


SZA 

♦ t-2 

8830 

ISOURCE 


JMP 

St art 1 

8890 

ISOURCE 


RET 

1 

8900 

ISOURCE 


END 

Dat a_trans 


8910 

8920 

8930 

8940 

8930 

8980 

8970 

8980 

8990 

9000 


' Get parameters of diti storage array. 
• Get parameters if data input array. 

1 Get number of samples. 


i 

i ***** 
i ** 
i ***** 
i 


Hist: 

ON KEY 
Nn*9 
Ns temp 


***************************************************************** * ' 
HISTOGRAM SUBPROGRAM ■> < 

******************************************************************* 

i Subroutine to produce online 
• histograms of raw LDV data. 

#0 GOTO 9270 


“Ns 


38 


4 


I 


%* 


a» 


9010 * 

9020 INPUT “No. of d it * samples per point <2000 maxm . > :",N» 

9O?0 ' 

9040 REDIN Data<3,Ns > 

9050 REDIN Dl<Nt,6> 

9060 DISP "Pr#i» CONT to initt*tt 6- ta Acquisition, press K0 to return to mein 


program. “ 

9070 PAUSE 
9080 GCLEAR 
9090 GRAPHICS 
9100 ! 

9110 RESET 10 

9120 CONTROL NASK 10; 1 

9130 WRITE 10 10.5J0 

9140 URITE 10 18,5} 1 1 stArt handshake by setting CTL0 

9150 Nt-6*Ns 

9160 FOR I '1 TO 5 

9170 Dummy-READBIN:t0> 

9180 NEXT I 

9190 ! 

9200 ENTER 10 WFHS Nt NOFORHAT; D 1 < *> ' f ASt dAt a ‘acqul s i 1 1 on 

9210 URITE 10 10,3;0 

9220 ICALL Daia trAnj 

9230 ICALL Dr Aw~h1 st 

9240 IF Ns«2000 THEN GOTO 9270 

9250 Nn-1 

9260 GOTO 9110 

9270 IF Ns-2000 THEN DUNP GRAPHICS 

9280 URITE 10 10.5J0 

9290 Ns m Ns temp 

9300 RED I N~Dat A< 3, Ns > 

9310 REDIN D 1 ( Ns , 6 ) 

9320 EXIT GRAPHICS 
9330 RETURN 


9340 

9350 

9360 

9370 

9390 

9390 

9400 

9410 

9420 

9430 

9440 

9450 


i 

ISOURCE 

ISCURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 


NAN Draw_hist 

The subroutine works by f'rst going to the sampl es-per-count 
tablet And using them to draw the old histograms in black 
<Bit**t)) to erase them. Then it calculates new samples-per- 
count tables from the data acquired from the LDV, and uses 
the new tables to draw histograms in white <Bit*l>. Then the 
subroutine returns to the main program. The very first time 
determined by Ntimes<“Nn>) the subroutine is called it dotin' 
erase th-? old histograms because there aren't any. 


9460 

ISOURCE 


EXT 

Get_val ue 

1 Declare 

subroutines stored 

9470 

ISOURCE 


EXT 

Get_i nf o 

1 outside 

of the main program. 

9480 

ISOURCE 


EXT 

Get_e 1 ement 



9490 

ISOURCE 

i 





9500 

ISOURCE 


CON 




9510 

ISOURCE 

Dat a_par! 

INT 

<*> 

1 Declare 

common variables. 

9520 

ISOURCE 

D l_p ar ; 

INT 

<*> 



9530 

ISOURCE 

Ns_par : 

INT 




9540 

ISOURCE 

Nt i me_p ir ; 

INT 




9550 

ISOURCE 

i 





9560 

ISOURCE 

Arrayd! 

BSS 

39 

1 Reserve 

spac » for data array 

9570 

ISOURCE 

El ementd: 

EQU 

Arr ayd* 1 6 

1 descriptor. 

9580 

ISOURCE 

Array 1 : 

ESS 

1024 

1 Reserve 

space for tables which hoto 

9590 

ISOURCE 

Array2: 

BSS 

1824 

■ the number of samples per count. 

9600 

ISOURCE 

Array3: 

BSS 

1024 




39 


9610 

ISOURCE 

Ns: 

BSS 

1 


9620 

ISOURCE 

Count : 

BSS 

1 1 Count 

and I are general purpose index vari- 

9630 

ISOURCE 

i: 

BSS 

1 'able* 

. Count is usually 0, 1, or 2, to denote 

9640 

ISOURCE 

Check : 

BSS 

1 i whether U.V, or U is being calculated. 

9650 

ISOURCE 

Int : 

BSS 

1 ' Int, 

Address, and Offset are all general 

9660 

ISOURCE 

Rddress : 

BSS 

1 I purpose storage areas. 

9670 

ISOURCE 

X_coord : 

BSS 

1 


9680 

ISOURCE 

Y coord: 

BSS 

1 


9690 

ISOURCE 

Bit: 

BSS 

1 


9700 

ISOURCE 

Nt 1 ires: 

BSS 

1 


9710 

ISOURCE 

i 




9720 

ISOURCE 


LIT 

30 


9730 

ISOURCE 

1 




9740 

ISOURCE 


SUB 



9750 

ISOURCE 

Draw_h i it : 

LOR 

-Ms 

• Get number of samples. 

9760 

ISOURCE 


LOB 

■Ns_par 


9770 

ISOURCE 


JSM 

Get_u**l ue 


9780 

ISOURCl 


LDH 

-Rrrayd 

1 Get parameters of data array. 

9790 

ISOURC' 


LDIl 

-Oat a_par 


9800 

ISOURCt* 


JSM 

Get_i nf o 


9810 

ISOURCE 


LDA 

-Nt i mes 

< Get Ntimes, which tells if this is 

9820 

ISOURCE 


LOB 

■Nt i me__par 

• the first time Draw_hist is being 

9830 

ISOURCE 


JSM 

Get_val ue 

1 cal 1 ed. 

9840 

ISOURCE 

i 




9850 

ISOURCE 


LOR 

-0 

• If this is the f.rst time Praw_hist 

9860 

ISCURCE 


STR 

Bit 

1 is being called, then jump to the 

9870 

ISOURCE 


LOR 

Nt i mes 

• data acquisition section. If not then 

9880 

ISOURCE 


RZR 

*+2 

1 write over the old histogram witn 

9890 

ISOURCE 


JMP 

Require 

1 black to erase it 

990O 

ISOURCE 

i 




9910 

ISOURCE 

flo_graph: 

LOR 

-0 

i Produce histograms for all three 

9920 

ISOURCE 


STR 

Count 

' channels. 

9930 

ISO'JRCE 

Mak e_h i it : 

LOR 

Count 

i Go to the end of the appropriate 

9940 

ISOURCE 


ROR 

-1 

' samples-pei — count table. 

9950 

ISOURCE 


LOB 

-1024 


9960 

ISOURCE 


STB 

I 


_aaza_ 

-ISOURCE 


MPY 



9980 

ISOURCE 


ROR 

-Array 1 


9990 

ISOURCE 


ROR 

— 1 


10000 

ISOURCE 


STR 

Address 


10010 

ISOURCE 

Mak e_rod: 

LOB 

Address , I 

1 Add together each adjacent pair of 

10020 

ISOURCE 


DSZ 

Address 

1 entries in the table. 

10030 

ISOURCE 


ROB 

Address, 1 


10040 

ISOURCE 


OSZ 

Address 


10050 

ISOURCE 


DSZ 

I 


10060 

ISOURCE 


RZB 

* + 2 

i If the result is zero, go to the next 

10070 

ISOURCE 


JMP 

Sk i p_b 1 1 s 

• pair of entries. If not, make sure 

10080 

ISOURCE 


STB 

Int 

i the sum is <150 and then draw a 

10090 

ISOURCE 


ROB 

— 150 

1 column with height equal to the sum. 

10100 

ISOURCE 


SBM 

i> + 3 


10110 

ISOURCE 


LOB 

"150 


10120 

ISOURCE 


STB 

Int 


10130 

ISOURCE 

i 




10140 

ISOURCE 


LDfl 

r 

• Calculate the X-cooroinate of the 

10150 

I SOURCE 


SRR 

1 

• column. 

10160 

ISOURCE 


ROR 

-20 

• 

10170 

ISOURCE 


STR 

.”_c oord 


10180 

ISOURCF 

i 




10190 

ISOURCE 

Mak e_b i : 

LOR 

"150 

* Calculate the Y-coordnate of the 

10200 

ISOURCE 


LDB 

Count 

1 top of the column. 
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<■ \ , 

S*!-J t *» 


10210 

ISOURCE 


ADB 

-1 

\ 


10220 

ISOURCE 


MPY 


kt V 


10230 

ISOURCE 


LOB 

Int 

k h 


10240 

ISOURCE 


TCB 




10230 

ISOURCE 


ADR 

B 

[■ * '' 


10260 

ISOURCE 


STR 

Y_coord 

: t- J 


10270 

ISOURCE 

1 



b'' V 


10280 

ISOURCE 


LOR 

-13 



10290 

ISOURCE 


STR 

Pa 

£* i 


10300 

ISOURCE 


LOR 

-31B 

if 


10310 

ISOURCE 


SFC 

* 


10320 

ISOURCE 


STR 

R3 

.5- ^ 


10330 

ISOURCE 


LOR 

Y ecord 

t ♦ 


10340 

ISOURCE 


LOB 

-36 


10350 

ISOURCE 


MPY 



10360 

ISOURCE 


LOB 

X coord 

i ? 


10370 

ISOURCE 


SBR 

4~ 

* ? 

* 

10380 

ISOURCE 

• 

RDR 

B 

f > 

A. ; 


10390 

ISOURCE 


CMR 


1 


10400 

ISOURCE 


SFC 

* 

? - 


10410 

ISOURCE 


STR 

R4 

fc- c 


10420 

ISOURCE 


STR 

R7 

* L 


10430 

ISOURCE 

1 



M* 
*» ? 


10440 

ISOURCE 


LOR 

X coord 


10430 

ISOURCE 


AND 

-I7B 

C~ $ 

i 


10460 

10470 

ISOURCE 

ISOURCE 


LOB 

SBL 

Bit 

15 



10480 

ISOURCE 


IOR 

B 



10490 

ISOURCE 


SFC 

* 


10500 

ISOURCE 


STR 

R4 

£ f 


10310 

ISOURCE 


STR 

R7 

f - 


10320 

ISOURCE 


DSZ 

Int 

% 


10530 

ISOURCE 


JMP 

Make_blt 

% ' 


10540 

ISOURCE 

1 



i 


10330 

ISOURCE 

Sk 1 p_bi - s: 

DSZ 

I 



10360 

ISOURCE 


JMP 

Make_rod 

I 


10570 

ISOURCE 

1 



| 


10580 

ISOURCE 


ISZ 

Count 



10390 

ISOURCE 


LOR 

—3 

r • 


10600 

isourcE 


RDR 

Count 

r ' 


10610 

ISOURCE 


SZR 

*+2 

i 


10620 

ISOURCE 


JMP 

Make_hi s 

• 


10630 

ISOURCE 

1 



L • 


10640 

ISOURCE 

flcqui rc : 

LOR 

Bit 

[ t 


10630 

ISOURCE 


SZR 

*+2 

r ' 


10660 

ISOURCE 


JMP 

Stop 



10670 

ISOURCE 

1 



S : 


10680 

ISOURCE 


LOR 

-flrrayl 

fc -* 


10690 

ISOURCE 


LOB 

-192 



10700 

ISOURCE 

Cl 

CLR 

16 

i 


10710 

ISOURCE 


RDR 

-16 



10720 

ISOURCE 


DSZ 

B 

l ? 


10730 

ISOURCE 


JMP 

Clear 

.a 


10740 

ISOURCE 

1 



* 

?* r 


10750 

ISOURCE 


LDfl 

-0 

sT ^ 


10760 

ISOURCE 


STR 

Count 

* ^ 


10770 

ISOURCE 


LDfl 

-5 

E~ « 


10780 

ISOURCE 


STR 

Check 

k f 


10790 

ISOURCE 


LDfl 

N* 



10800 

ISOURCE 


STR 

I 
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Prepare the graphics screen for 
< nput . 


Calculate word location. 


Input word location. 
Calculate bit location. 

Input bit location. 


Make histograms for all three 
channel s. 


If Bit-1 then the new histograms 
have been drawn. Go back to the 
main program. If not, continue. 

Zero out the tables which hold the 
number of samples per count. 


Prepare to write new samples-per- 
count tables. 


* 




10810 

I SOURCE 

1 



10820 

ISOURCE 

Begin: 

LDA 

Count 

10830 

IS0URCE 


LDB 

Ns 

10840 

IS0URCE 


MPY 


10830 

ISOURCE 


ADA 

I 

13880 

IS0URCE 


ADA 

»- 1 

10870 

IS0URCE 


STA 

El emcntd 

10880 

IS0URCE 


LDA 

-Int 

10890 

I SOURCE 


LDB 

-Arrayd 

10908 

IS0URCE 


JSM 

Get_el enent 

1091 * 

I SOURCE 


LDA 

Int 

10920 

IS0URCE 


LDB 

-1023 

10930 

IS0URCE 


AND 

B 

10940 

IS0URCE 


LDB 

Count 

10930 

IS0URCE 


TCB 


10980 

IS0URCE 


ADB 

Check 

10970 

IS0URCE 


SZB 

Strai ght 

10980 

IS0URCE 


TCA 


10990 

IS0URCE 


ADA 

■1324 

11000 

IS0URCE 

Strai gh-, : 

STA 

Int 

11010 

IS0URCE 


LDA 

-1024 

11020 

I SOURCE 


LDB 

Count 

11030 

IS0URCE 


MPY 


11040 

IS0URCE 


ADA 

-Array l 

11050 

IS0URCE 


ADA 

Int 

11060 

IS0URCE 


ISZ 

A, I 

11070 

IS0URCE 

i 



11080 

IS0URCE 


ISZ 

Count 

1 1090 

IS0URCE 


LDA 

—3 

11100 

IS0URCE 


ADA 

Count 

11110 

I SOURCE 


SZA 

*+2 

11120 

IS0URCE 


JMP 

Begi n 

11130 

IS0URCE 


STA 

Count 

11140 

IS0URCE 

i 



1 1 150 

ISOURCE 


DSZ 

I 

11160 

ISOURCE 


JMP 

Beg l n 

11170 

ISCURCE 

i 



11180 

ISOURCE 


ISZ 

Bit 

11190 

ISOURCE 


JMP 

Do graph 

11200 

ISOURCE 

Stop: 

RET 

1 

11218 

ISOURCE 

• 



11220 

ISOURCE 


END 

Draw hist 


< Figure out which element of the 
• data array we want to pick up. 


1 Cet a raw datum from the data array. 


• Strip off the first six 

1 b'ts of the raw data word. 

* See if Count * Check. 


' If true, use the modified data 
< word as an index. If not, use 
' 1024 minus the data word, 

i Store the count we have gotten. 
' Increment the appropriate place 
• in the samp I es-per-c ount table 
1 by one. 


• Have U,V, and W all been done'’ 

■ If not, go back again. 

1 If so, set Count-0. 

1 Continue to fill sampl es-per-c ount 
1 tables until out of samples. 

< Now set Bit-1 and go back and d'au 
1 the new histograms in white. 



/ 


/ 
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REM PROGRAM STAT 

I PROGRAM TO REDUCE RAW DATA FROM THREE-COMPONENT LDV SYSTEM 
I Input is 3 channels of raw LV data and calibration constants 
I Output Is the three components of mean velocity) all the 
• components o- ^•'•'olds stress and some selected third-order 
I products. Can * --r out counts which are excessively far 

1 from the mean. display histograms of raw and 

! filtered dat i f«. .11 3 channels. Plots data at end and displays 

I reduced data for entire run. Can write a summary file containing 
! all reduced data and calibration constants. 

OPTION BASE 1 
ICOM 18000 
IDELETE ALU 
IRSSEMBLE Find vel 

COM INTEGER Dat - »< 3 , 2000> , N* , Ngs , Count b l n< 3, 1 024 ) 

COM INTEGER Rangel , Range2, Range3 

COM REAL Sdev, D* 1 , Df 2, Df 3 

COM REAL Theta f Nub, Numi xl , Numi x2, Numi x3 

COM REAL Su, Sv, 3w, Suu, Svv, Sww, Suv, Suu, Svw 

COM REAL Suuv, Suvv, Suuw, Suww, Suvw 

I 

DIM LK3>,R1 <3> , L 1 1<3),U1 1 <3> , Seal e<3> 

DIM DatefC 801 , F i 1 e*C63 , Name*C43 , Ti 1 1 #C303 

! 

REAL Xpos, YpoS) 2pos 
REAL Poi nt <20> 15> 

REAL Re, Ue 


INTEGER Fi Ino,A,Run 

f 

DATA 0, 0,0,0 

READ Max_y, Mi n_v, FI 1 ter, Sdev 
DATA 75.64,43.64, 15.64,99. ,69. ,39. 

MAT READ LI 1 , U 1 l 
DATA 5,46,87,41 .82, 123 
MAT READ Ll ,R1 

DATA -. 8, 2. 4, 0. , .033, -.0023, .01 

READ Llim<l>,Ulim<l>,Llim<2>,Ulim<2>,Llim<3>,Uf*m<3> 
DATA .4, .003,. 8023 
MAT READ Scale 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


" ** << PROGRAM STAT I 3-COMPONENT VELOCITY DATA >> ** 
“PROGRAM STRUCTURE" 

"1. Read raw data from floppy disc" 

"2. Convert to velocity" 

“3. Calculate statistics" 

"4. Print results" 

“5. Write to disc f i le" 


' ** Read raw data from floppy disc ** 

l 

INPUT “Enter pa"ent filename (or E to exit program) :“,Name$ 

IF Namef»“E“ THEN GOTO 1780 

INPUT “Enter no of first and last data files “,Filel,Nf 
Ki 1 l_u-3 

INPUT “Is this a twc-channel or three-channel run ( 2 ' 3 , Default 3> 


Hi st*-“Y“ 


T 


K i I 
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610 INPUT "Do you wish to see histograms of the data points <Y^N, Default Y> 
", H( st * 

620 Filter*-“N" 

630 INPUT "Do you wish to filter th* output <Y/N, Default N) Filter* 

640 IF Filtrl-'N' THEN GOTO 660 

650 INPUT "En-,er standard deviation to throw out",Sdev 

660 DISP “READING RAW DATA FROM DISK" 

670 MASS STORAGE IS ":H8,0,0" 

680 PRINTER IS 0 

690 F i I no»F il el 

700 IF F i 1 no>Nf THEN GOTO 1740 

710 Fi leS«Name*3,VAL.KFi lno> 

720 ASSIGN Fi le* TO #1 
730 ' 

740 READ #l?Date* 

750 READ # 1 ; T 1 1 1 * 

760 READ #i;Name* 

770 READ *l*,Dn 

780 READ #1 J Nub, Num i xl , Numi x2, Numi x3, That a, Run 
790 READ #i;Vref,Ue,Dfl,Df2,Df3 

800 READ *1 ; Xpos, Ypos, Zpos, Ns ' 

810 MAT READ #1 ;Dat i 
820 ' 

830 PRINT Tltl* 

840 PRINT "DATE AND TIME OF TEST Dates 
850 PRINT “FILE NAME ON DISK I", Name* 

860 PRINT "DISK NUM3ER : “ , Dn 

870 PRINT "RUN NO. :P0INT NO. IN PROFILE) i " , Run 
880 PRINT "X,Y,Z * ' , Xpos , Ypos, Zpos 
890 • 

900 ' CALCULATE VELOCITIES 

910 DISP “CALCULATING VELOCITIES AND TAKING RUNNING SUMS" 

920 ' 

930 ICALL Find vel 
940 • 

930 ' CALCULATE STATISTICS 

960 DISP "CALCULATING STATISTICS" 

970 ' 

989 IF Kill u*3 THEN GOTO 1030 

990 DATA 0,0, 0,0, 0,8,0 
1000 RESTORE 990 

1010 READ Sw, Suw , Suw , Sww , Suuw, Suww, Suvw 
1020 ' 

1030 Ubar^Su'Ngs 
1040 Vbar»Sw^Ngs 
1030 Wbar^Sw'Ngs 
1060 Upr i 2*Suu''Ngs-Ubar*Ubar 
1070 Vpr l 2»Svv'Ngs-Vbar*Vbar 
1080 Wpr i 2*Sww«'Ngs-Wbar'*Wbar 
1090 Uubar“Suw''Ngs-Ubar*Vbar 
1100 Uwbar «Suw/Ngs-Ubar*Wbar 
1110 Vwbar«Suw''Ngs-Vbar*Wbar 

1120 Uuvbar»$uuv/ / Ngs-2*Ubar *Suv/'Ngs-Vbar *Suu/Ngs+2*Vbar*Ubar*Ubar 

1130 Uuubar»Suuu^Ngs-2*Vbar*Suu^Ngs-Ubar-*Svv^Ngs+2*Ubar'»Vbar*Vbar 

1140 Uuubar»Suuw.- , Ng3-2*Ubar*Suw/'Ngs-Wbar*Suu/Ng3 + 2*Wbar'*Ubar*Ubar 

1150 Uwwbar»Suww^Ngs-2*Ubar*Suw^Ngs-Ubar*Sww/Ngs+2*Ubar *Wbar*Wbar 

1160 Uvwbar *Suuw^Ngs-Wbar ♦Suu/Ngs-Ubar *Sww^ Ngs- Vbar*Suw/Ngs+2*Ubar * Vbar *Wbar 

1170 ' 

1130 ' STORE RESULTS IN ARRAY FOR PLOTTING 

1 190 ' 

1200 F l*Fi 1 no-Fi I e 1 + l 


* 
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I 


l 



* 


1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1 otr* 
1770 
1780 
1790 
1800 


Point(Fl, 1 >*Ypos 

Point (FI ,2>*Ubar 

Point (FI ,3>*Vbar 

Pol nt (FI , 4 ) *Wbar 

Pot nt <F1 , 5>*Upr i 2 

Po i nt (FI , 6)**Vpr i 2 

Pol nt (FI , 7>*Upr i 2 

Point ( F 1 , 8 ) *Uvb tr 

Pol nt <F1 , 9)"Uub ir 

Point(Fl,10)-Vwbar 

Po i nt ( F 1 , 1 1 > »Uuvbar 

Poi nt <F1 , 12)»Uvvbar 

Pomt(Fl, 13)»Uuwbar 

Poi nt <F1 , 14>»Uwwbar 

Poi nt <F1 , 15>»Uvubar 

IF Ypos<Max_y THEN GOTO 1390 

Max_y*Ypos 

Min u«Ubar 

IF Ypo»>Mi n_y THEN GOTO 1450 

Mi n_y»Ypos 

Max_u«Ubar 

! 

' PRINT RESULTS 

i 

PRINT 

PRINT "Bragg shift ft equency * “ , Nub 
PRINT “Mixing f requenc l es B " , Numi x 1 , Num l x2 , Num i x3 
PRINT “Third beam angle ■ ", Theta 
PRINT “Run number ■ “,Run 

PRINT "Tunnel reference voltage (volts) * “,Vref 

PRINT “Fringe spacings ■ " , Df 1 , Df 2 , Df 3 

PRINT "No. of samples = “,Ns 

PRINT "Ranges ” ", Rangel , Range2, Range3 

PRINT “Filter Clio, of std. deviations! * “,Sdev 

PRINT "Adjusted no. of samples *“,Ngs 

PRINT "Ubar - “ , *R0UND(Ubar,-4> 

PRINT “‘oar ■ " ,PROUNI»(Vbar, -4) 

PRINT "Ubar - “ 1 PR0UNB(Wbar,-4-i 
PRINT “Upr i 2 - ' , PROUND(Upr l 2 , -4 ) 

PRINT “Vpr i 2 - * , PROUND ( Vpri 2 , -4 > 

PRINT “Wpr l 2 - * , PROUND (Wpr t 2, -4) 

PRINT “ UVbar - ' , PROUND(Uvbar, -4) 

PRINT "UWbar « ' , PROUND(Uwbar , -4) 

PRINT " VWbar ■ ',FR0UND(Vubar,-4> 

PR T NT “UUVbar * “ , PROUND(Uuvbar, -5) 

PRINT "UWbar - ‘ , PROUND(Uvvbar , -5) 

PRINT "UUWbar - " , PROUND(Uuubar , -5) 

PRINT "UUWbar « " , PR0UND( Uwwbar , -5 > 

PRI‘!T “UVWbar - “ , PROUND(Uvwbar , -5) 

PRINT 

IF HistS-"Y" THEN GOSUB Histogram , 

Fi lno-Fi lno+l 
GOTO 700 
Ue*Max_u-M l n_u 
Plotr*-“Y" 

INPJT " Do you wish to see the turbulence quantities 

IF P1otr»«"Y" THEN GOSUB Plot 
MASS STORAGE IS “:H8,0,1" 

PRINTER IS 16 
END 


(Y/N, 


Default Y > '>",P 
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1310 ! 

1820 ' *•«**« *******-«#*******P1.0T DATA****************4************************ 

1830 ' 

1840 Plot: PRINTER I *3 16 

1850 PRINT "Estimate of normalizing velocity .s :";Ue 

1860 INPUT "Enter correct normalizing velocity if di fferent . " , Ue 

1370 Graph*»“Y“ 

1880 INPUT “Do you wish to graph the turbulence quantities <YxN, Default Y) ">", 
Graph* 

1830 PRINTER IS 0 

1900 PRINT "Normalizing velocity is :";Ue 
1910 ! 

1920 FOR 1-1 TO 4 
1930 DATA l, 2, 2, 3 
1940 READ Exp 
1950 Kn-3 

1960 IF I *4 THEN Kn*5 
1970 FOR J«1 TO Nf-N*' i 1 e+1 
1980 FOR K-l TO Kn 
1990 Index»d-1 )*3+K+l 

2000 Polnt<J, Index) *Poi nt < J, Index)/Ue , 'Exp 

2010 IF Index>10 THEN Poi nt < J , Index ) »Po i nt < J , Index)*lO00 

2020 NEXT K 

2030 NEXT J 

2040 NEXT I 

2050 IF Graph*>"N" THEN GOTO 2860 
2060 • 

2070 PLOTTER IS 13, “GRAPHICS" 

2080 GRAPHICS 
2090 DEG 
2100 • 

2110 Ymax*Poi nt < 1 , 1 ) 

2120 Ym i n*Ymax 

2130 FOR J-l TO Nf-Filel+1 

2140 IF Pot nt < J, 1 ) >Ymax THEN Ymax*Poi nt < J , 1 ) 

2150 IF Point<J,l)<Yi#in THEN Ym i n*Po l nt < J , l > 

2160 NEXT J 
2170 ' 

2130 FOR 1=1 TO 3 

2190 LIMIT 0, 184.47,0, 149. » 

2200 LOCATE LI <I),R1 :i),39,99 
2210 CLIP LI <I),R1 <1 >,39,99 
2220 ' 

2230 Xmax*Point < 1 , < 1-1 )*3*2) 

2240 Xmin*Xmax 

2250 FOR J-l TO Nf-Filel+1 

2260 FOR K* 1 TO 3 

2270 IF Poi nt < J, < 1-1 >*3+K+l ) >Xmax THEN Xmax=Po i nt < J , < I - 1 ) *3+K+ 1 ) 

2280 IF Pol nt < J, <1-1 >*3+K+i)<Xmin THEN Xmtn=»Point(J, d-l)43+K + l) 

2290 NEXT K 
2300 NEXT J 
2310 ' 

2320 SCALE L 1 i m< I ) , U I i m< I ) , Ym i n, Ymax 
2330 FRAME 

2340 AXES Scaled), . 1,8,0 
2350 UNCLIP 
2360 ' 

2370 LORG 3 
2380 CSIZE 2.5 
2393 Ypos*Ymin 

2400 MOVE LI imd)+Sc t1ed)^5,Ypos 
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2410 DRAM LI i »< I ) , Ypos 
2420 SETGU 
2430 RPLOT 2,0, -2 
2440 SETUU 

2450 LABEL DR0UND< Ypos, 2) 

2460 Ypos*Ypos+. 1 

2470 IF Ypo»< "Y biax THEN GOTO 2400 
2480 ! 

2490 LORG 6 
2500 Xpos=Ll > i»< I ) 

2510 MOVE Xpos, Ymt n- . 04 
2520 LABEL DROUND<Xpos, 2> 

2530 MOVE Xpos, Ymi n 
2540 DRAM Xpos, Ymi n+ , 05 
2550 Xpos*Xpos+Sc*l * : I ) 

2560 IF Xpos<Ul i m< I ? THEN GOTO 2510 
2570 ! 

2580 FOR J=1 TO Nf-Fi 1*1+1 
3£90 Ypos*Poi nt < J, 1 > 

2600 FOR K=1 TO 3 

2610 MOVE P01nt<J,<I-l)*3+K+l),Ypos 
2620 SETGU 

2630 ON K GOTO 2640, 2700, 2760 

2640 FOR Arc»0 TO 360 STEP 20 

2650 PDIR Arc 

2660 RPLOT .5,0 

2670 NEXT Arc 

2680 PDIR 0 

2690 GOTO 2810 

2700 RPLOT .5, .5,-2 

2710 RPLOT .5, -.5,-1 

2720 RPLOT -.5, -.5,-1 

2730 RPLOT -.5, .5,-1 

2740 RPLOT .5, .5,-1 

2750 GOTO 2810 

2760 RPLOT 0, .5,-2 

2770 RPLOT 0, -. 5, -1 

2780 RPLOT . 5, 0, -2 

2790 RPLOT -.5, 0,-1 

2800 GOTO 28 1 C 

2810 SETUU 

2820 NEXT K 

2830 NEXT J 

2840 NEXT I 

2850 DUMP GRAPHICS 

2860 IMAGE 2X"Y"10X, *u2-'Uo2“ , 3X, , •v2✓Uo2 ,, , 3X, "u2/Uo2" , 5X, “uv^Uo2" , 3X, "uu-'Uo2" , 3X 
, "vw/Uo2",/ 

2870 PRINT USING 2860 
2880 FOR 1=1 TO Nf-F i 1*1+1 

2890 IMAGE MDD. DD , 6X , MD. DDDD ; 2X , MD . DDDD , 2X , MD . DDDD , 4X , MD . DDDD , 2X, MD . D.DDD , 2X , MD . 
DDDD 

2900 PRINT USING 2898; Point <I,l),Point<I,5>, Point (I, 6>,Point(I, Point (I, 8>,Po 
1 nt < 1 , 9 > , Po l nt < I , 10) 

2910 NEXT I 

2920 IMAGE 2',3X,"Y" , 9*., "uciv/Ij 2 « , JX, -wvv/Uo3" , 5X, "uuw/Uo3 M , 3X, " uww/ Uo3" , 5X, M u 
M(j/Uo3" , / 

2930 PR T NT USING 2928 
2940 FOR 1=1 TO Nf -F i 1 *♦ 1 

2950 IMAGE MDD. DD, 6X , MD. 5D, 2X, MD. 5D, 4X, MD. 5D, 2X, MD. 5D, 4X, MD. 5D 

2960 -'MNT USING 2958; Pc i nt < 1 , 1 ) , Pot nt < 1 , 1 1 ) , Po i nt < 1 , 1 2 ) , Po 1 nt < 1 , 1 3 > , Po \ nt < 1 , 1 4 
),Poi. 1,15) 



t 


I 


2970 NEXT I 
2980 i 

2990 ' WRITE SUMMARY DATA FILE 

3000 I 

3010 Sum*-"N“ 

3020 INPUT “Bo you wish to writ* a Summary Data File <Y/N Default N) ^"fSumS 
3030 IF SumS- - N” THEN GOTO 3210 
3040 Filett-Name* 

3030 DISP "File "jFile*;" being written to disk" 

3060 MASS STORAGE IS ":H8,0,0“ 

3070 CREATE Fi le*,40 
3080 ASSIGN Fi 1e* TO #1 
3090 PRINT ttl ; Date* 

3100 PRINT tttjTitl* 

3110 PRINT ttl ; Name* 

3120 PRINT ttl ; Dn 

3130 PRINT ttl ; Nub, Num i xl , Num i x2, Numt x3 , Thet a, Run 

3140 PRINT ttl ; Vref , U*,Dfl,Df2,Df3 

3130 PRINT ttlJNs.Fi Wl,Nf 

3160 FOR 1*1 TO Nf-Nfile+1 

3170 FOR J-l TO 13 

3180 PRINT *1 J Pol nt < t , J> 

3190 NEXT J 
3200 NEXT I 
3210 RETURN 
3220 I 

3230 1 *********************DRAW HISTOGRAMS***********#***-********************* 

3240 ! 

3230 Histogram: PLOTTER IS 13, "GRAPHICS- 

3260 GRAPHICS 
3270 K-l 

3280 FOR 1*1 TO 3 
3290 Cmax-0 

3300 FOR J-l TO 1023 STEP 2 

3310 Countbm<I,J)-Countbin<I, J> +Count bin<I,J+l> 

3320 IF Countbi n< I , J > >Cmax THEN Cmax-Countbi n( I , J> 

3330 NEXT J 

33 °) IF <Cmax-0> OR :Cmax-2000> THEN GOTO 3660 

3350 LIMIT 0, 184.47,13, 149.8 

3360 LOCATE 7. 03, 1 19 . 63, LI 1 <IO , J1 1 <K> 

3370 CLIP 7.03, 119. 63, LI l(K),UU(IO 
3380 K*K*l 

3390 SCALE l , 1024, 0, Cmax 

3400 FRAME 

3410 UNCLIP 

3420 LINE TYPE 1 

3430 CSIZE 2.8 

3440 LORG 8 

3430 FOR Y-l TO 4 

3460 Ypos*Y*Cmax/'4 

3470 PLOT 10, Ypos, -2 

3480 PLOT -9, Ypos, -1 

3490 MOVE 3, Ypos 

3500 LABEL PROUND < Ypos , to > 

3510 NEXT Y 

3520 LORG 6 

3330 Ypos*Ciiax''20 

3340 FOR X*i TO 8 

3350 Xpos-X* 1 28 

3360 PLOT Xpos, Ypos, -2 
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3370 

PLOT Xpos, -Ypos , -1 


3588 

MOVE Xpos, -Ypos*l . 

3 

3590 

LABEL Xpos 


3600 

NEXT X 



3610 

FOR J«1 

TO 1023 STEP 2 

3620 

IF Countbi n C I , J >«0 

THEN GOTO 3650 

3630 

MOVE J, 

9 


3640 

DRAM J, Countbi n : I , 


3650 

NEXT J 



3660 

NEXT I 



3670 

MOVE 512, -Ypos*5 


3680 

LORG 6 



3690 

LABEL “ 

PLOT OF COUNT VS NUMBER OF SAMPLES PER COUNT" 

3700 

DUMP GRAPHICS -Ypos*10 

3710 

RETURN 



3720 

! 



3730 

1 



3740 

ISOURCE 


NAM Find_vel 

3750 

ISOURCE 

1 


3760 

ISOURCE 

1 This jubroutine converts raw data counts to 1 nstantaneous 

3770 

ISOURCE 

• velocities, then sums several different products of the 

3780 

ISOURCE 

1 velocity 

components. All i nput and output data is passed 

3790 

ISOURCE 

1 through 

the COMMON stoarage area. The inputs are the raw 

3800 

ISOURCE 

• data trray (Arrayd), the Bragg shift frequency <Nub>, the 

3810 

ISOURCE 

1 mixing frequencies (Numixl, Numix2, Numix3>, the fringe 

3820 

ISOURCE 

' spaemgs 

<Dfl, Df2, Df3>, the crossing angle of the third 

3830 

ISOURCE 

1 beam CTheta), and the number of samples in a data point (N 

3840 

ISOURCE 

' The outputs are the summations of various products of the 

3850 

ISOURCE 

1 veloc i ty 

components, including U, V, W, U*U, V*V, W*U, 

3860 

ISOURCE 

' U*V, U*W 

, V*W, 0*U*V, U*V*V , U*U*W, U*W*W, and U*V*W. 

3870 

ISOURCE 

i 


3880 

ISOURCE 


EXT Get_value 1 Declare subroutines stored 

3890 

ISOURCE 


EXT Get_info ' outside of the main program. 

3980 

ISOURCE 


EXT Get_e lament 

3910 

ISOURCE 


EXT Put_element 

3920 

ISOURCE 


EXT Int_to_rel 

3938 

ISOURCE 


EXT Rel_to_int 

3940 

ISOURCE 


EXT Rel_math 

3950 

ISOURCE 


EXT Put_val ue 

3960 

I SOURCE 

i 


3970 

ISOURCE 


COM 

3980 

ISOURCE 

Data_par: 

INT <*> 1 Declare common variables. 

3990 

ISOURCE 

Ns_par : 

INT 

4000 

ISOURCE 

Ngs_par : 

INT 

4010 

ISOURCE 

Cb i n_par: 

INT <#> 

4020 

ISOURCE 

Ranl_par : 

INT 

4030 

ISOURCE 

Ran2_par : 

INT 

4040 

ISOURCE 

Ran3_par : 

INT 

4050 

ISOURCE 

Sdev_par : 

REL 

4063 

ISOURCE 

Dfl par: 

REL 

4070 

ISOURCE 

Df2_par: 

REL 

4080 

ISOURCE 

Df 3_par : 

REL 

4090 

ISOURCE 

Theta_p »r: 

REL 

4100 

ISOURCE 

Nub_par : 

REL 

4110 

ISOURCE 

Nm i x 1 _p »r : 

REL 

4120 

ISOURCE 

Nmi x2_p u : 

REL 

4138 

ISOURCE 

Nml x3_p tr : 

REL 

4140 

ISOURCE 

Su_par : 

REL 

4150 

ISOURCE 

Sv_par : 

REL 

4160 

ISOURCE 

Sw_par : 

REL 



\ 
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4170 

ISOURCE 

Suu_par : 

REL 



4180 

ISOURCE 

SwjliT : 

REL 



4190 

ISOURCE 

Swu_p ar : 

REL 



4200 

ISOURCE 

° iv_par : 

REL 



4210 

ISOURCE 

Suu_par : 

REL 



4220 

ISOURCE 

Svu_par : 

REL 



4230 

ISOURCE 

Suuv_par : 

REL 



4240 

ISOURCE 

Suvv_par : 

REL 



4250 

ISOURCE 

Suuu_par : 

REL 



4260 

ISOURCE 

Suwu_par : 

REL 



4270 

ISOURCE 

Suvu ptri 

REL 



4280 

ISOURCE 

i 




4290 

ISOURCE 

Arrayd: 

BSS 

39 

1 Reserue space for data array 

4300 

ISOURCE 

El ementd: 

EQU 

Arrayd+16 1 descriptor. 

4310 

ISOURCE 

Array 1 : 

BSS 

4096 

1 Reserue space for lookup tables 

4320 

ISOURCE 

Array2: 

BSS 

4096 

1 for count to uelocity couersion. 

4330 

ISOURCE 

flrray3: 

BSS 

4096 


4340 

ISOURCE 

B l n_u : 

BSS 

1024 

1 Reserue space '‘or table of counts 

4350 

ISOURCE 

B l n_v 1 

BSS 

1024* 

1 used to generate histograms. 

4360 

ISOURCE 

Bi n_u: 

BSS 

1024 


4370 

ISOURCE 

Rangel : 

BSS 

1 

1 Reserue space for uarlous Input 

4380 

ISOURCE 

Range2: 

BSS 

1 

1 and output uariables. 

4390 

ISOURCE 

Range3: 

BSS 

1 


4400 

ISOURCE 

Ns: 

BSS 

1 


4410 

ISOURCE 

Df l: 

BSS 

4 


4420 

ISOURCE 

Df2: 

BSS 

4 


4430 

ISOURCE 

Df3: 

BSS 

4 


4440 

ISOURCE 

Theta: 

BSS 

4 


4450 

ISOURCE 

Nub: 

BSS 

4 


4460 

ISOURCE 

Numi xl : 

BSS 

4 


4470 

ISOURCE 

Nuiji x2: 

BSS 

4 


4480 

ISOURCE 

Numi x3: 

BSS 

4 


4490 

ISOURCE 

Su: 

BSS 

4 


4500 

ISOURCE 

Su: 

BSS 

4 


4519 

ISOURCE 

Su: 

BSS 

4 


4520 

ISOURCE 

Suu: 

BSS 

4 


4530 

ISOURCE 

Suu: 

BSS 

4 


4540 

ISOURCE 

Suu: 

BSS 

4 


4550 

ISOURCE 

Suu: 

BSS 

4 


4560 

ISOURCE 

Suu: 

BSS 

4 


4570 

ISOURCE 

Suu: 

BSS 

4 


4580 

ISOURCE 

Suuu: 

BSS 

4 


4590 

ISOURCE 

Suuu: 

BSS 

4 


4600 

ISOURCE 

Suuu: 

BSS 

4 


4610 

ISOURCE 

Suuu: 

BSS 

4 


4620 

ISOURCE 

Suuu: 

BSS 

<< 


4630 

ISOURCE 

Count : 

BSS 

i t 

Count and I are general purpose index uari- 

4640 

ISOURCE 

I: 

BSS 

1 i 

ables. Count is usually 0, 1, or 2 to denote 

4650 

ISOURCE 

Check : 

BSS 

1 ' 

whether U, V, or W -s being calculated. 

4660 

ISOURCE 

Int : 

BSS 

1 ' 

Int, Address, <’nd Offset are all general 

4670 

ISOURCE 

In-2: 

BSS 

1 ' 

purpose storage areas. 

4680 

ISOURCE 

Address : 

BSS 

1 


4690 

ISOURCE 

Offset : 

BSS 

l 


4780 

ISOURCE 

R 1 : 

BSS 

4 ' 

Rl, R2, and R3 are the count-to- 

4710 

ISOURCE 

R2: 

BSS 

4 i 

frequencv conversion factors. 

4720 

ISOURCE 

R3: 

BSS 

4 

• 

4730 

ISOURCE 

Xvar : 

BSS 

4 ' 

Xuar and Yvar are general purpose real 

4740 

ISOURCE 

Vuar : 

BSS 

4 ' 

number storage »rets. 

4750 

ISOURCE 

u: 

BSS 

4 ' 

U, V, and W are the instantaneous uelocity 

4760 

ISOURCE 

v: 

BSS 

4 ' 

c oriponent s . 


8 
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I 


X 


I 


4770 

ISOURCE 

w: BSS 

4 

4780 

I SOURCE 

Uu: bss 

4 

4790 

ISOURCE 

Uv: bss 

4 

4800 

ISOURCE 

Cos: BSS 

4 * Cos ■ 

<810 

ISOURCE 

sin: BSS 

4 

4820 

ISOURCE 

Max_u: BSS 

1 

4830 

ISOURCE 

n«x_u: B~ 3 

1 

4840 

ISOURCE 

Max_w: BSS 

1 

4850 

ISOURCE 

Avg_u: BSS 

4 

4860 

ISOURCE 

Hug_u: BSS 

4 

4870 

ISOURCE 

Avg u: BSS 

4 

4880 

ISOURCE 

Dev~u: BSS 

4 

4890 

ISOURCE 

Dev_v: BSS 

4 

4900 

ISOURCE 

Dev_w: BSS 

4 

4910 

ISOURCE 

Ngs: BSS 

1 

4920 

ISOURCE 

Sdev: BSS 

4 

4930 

ISOURCE 

Rad: DAT 

5. 729578E1 

4940 

ISOURCE 

Mill: DAT 

1.E6 

4950 

ISOURCE 

One: DAT 

1. 

4960 

ISOURCE 

Zero: DAT 

0. 

4970 

ISOURCE 

LIT 

80 

4980 

ISOURCE 

i 


4990 

ISOURCE 

SUB 


5000 

ISOURCE 

Find ve 1 : LDA 

■Ns 

5010 

ISOURCE 

LDB 

=Ns_par 

5020 

ISOURCE 

JSM 

Get_va! ue 

5030 

ISOURCE 

LDA 

■Arr ayd 

5040 

ISOURCE 

LDB 

■ Dat a_par 

5050 

ISOURCE 

JSM 

Get_i nf o 

5060 

ISOURCE 

LDA 

■Sdev 

5070 

ISOURCE 

LDB 

■Sdev_par 

5080 

ISOURCE 

JSM 

Get value 

5090 

ISOURCE 

LDA 

■Df 1 

5100 

ISOURCE 

LDB 

■Df l_par 

5110 

ISOURCE 

JSM 

Get value 

5120 

ISOURCE 

LDA 

■Df 2 

5138 

ISOURCE 

LDB 

■Df 2_par 

5140 

ISOURCE 

JSM 

Get value 

5150 

ISOURCE 

LDA 

■Df 3 

5160 

ISOURCE 

LDB 

■Df 3_par 

5170 

ISOURCE 

JSM 

Get_va 1 ue 

5180 

ISOURCE 

LDA 

■ Thet a 

5190 

ISOURCE 

LDB 

■ Thet a_par 

5200 

ISOURCE 

JSM 

Get_val ue 

5210 

ISOURCE 

LDA 

■Nub 

5220 

ISOURCE 

LDB 

■Nub_par 

5230 

ISOURCE 

JSM 

Get_val ue 

5240 

ISOURCE 

LDA 

■Numi x 1 

5250 

ISOURCE 

LDB 

■Nmi xl_par 

5260 

ISOURCE 

JSM 

Cet_val ue 

5270 

ISOURCE 

LD'i 

■Numt x2 

5280 

ISOURCE 

LDB 

■Nmi x2_par 

5290 

ISOURCE 

JSM 

Get_val ue 

5300 

ISOURCE 

LDA 

■Numi x3 

5310 

ISOURCE 

LDB 

■Nmi x3_par 

5320 

ISOURCE 

JSM 

Get_val ue 

5330 

ISOURCE 

l 


5340 

ISOURCE 

' ihe loop headed by Get_' 

5350 

ISOURCE 

1 the count-to- 

-frequency < 

5360 

ISOURCE 

• the range) for U,V, and 


Cos *nd Sin are the cos and sin of Theta. 


' Get number of samples. 


Get parameters of data a~ray. 


■ Get input parameters. 


q is repealed three times to get 
conversion factors <which depend on 
W. Whenever a loop is controlled by 
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5370 

ISOURCE 

i the variable 

"Count'*, th 

5380 

5390 

5400 

ISOURCE 

ISOURCE 

ISOURCE 

' the sune 

i 

i 

for 

U,V, and U. 

5410 

ISOURCE 


LDfl 

-0 

5420 

ISOURCE 


STA 

Count 

5430 

ISOURCE 

Get_freq: 

LDfl 

Count 

5440 

5450 

ISOURCE 

ISOURCE 


LDB 

MPY 

Ns 

5460 

ISOURCE 


STfl 

E 1 ement d 

5470 

ISOURCE 


LDfl 

■ Int 

5480 

ISOURCE 


LDB 

■flrrayd 

5490 

ISOURCE 


JSM 

Get_e 1 ement 

5500 

ISOURCE 


LDfl 

Int 

5510 

ISOURCE 


LDB 

-15360 

5520 

ISOURCE 


AND 

B 

5530 

5540 

ISOURCE 

ISOURCE 


SflR 

TCfl 

10 

5550 

ISOURCE 


LDB 

- 15 

5560 

ISOURCE 


ADR 

B 

5570 

ISOURCE 


LDB 

■Rangel 

5580 

ISOURCE 


ADB 

Count 

5590 

ISOURCE 


STfl 

B.I 

5600 

ISOURCE 


LDB 

-1 

5610 

ISOURCE 


SZfl 

Loopend 

5620 

ISOURCE 

Loop: 

SBL 

1 

5630 

ISOURCE 


DSZ 

A 

5640 

ISOURCE 


JMP 

Loop 

5650 

ISOURCE 

Loopend : 

STB 

Int 

5660 

ISOURCE 


LDfl 

-Int 

5670 

ISOURCE 


STfl 

Oper__l 

5680 

ISOURCE 


LDfl 

-Yvar 

5690 

ISOURCE 


STfl 

Resul t 

5700 

ISOURCE 


JSM 

Int to rel 

5710 

ISOURCE 


LDfl 

— 3.2E4 

5720 

ISOURCE 


LDB 

-Xvar 

5730 

ISOURCE 


XFR 

4 

5740 

ISOURCE 


STB 

Oper_l 

5750 

ISOURCE 


LDfl 

-Yvar 

5760 

ISOURCE 


STfl 

0per_2 

5770 

ISOURCE 


LDfl 

Count 

5780 

ISOURCE 


SflL 

2 

5790 

ISOURCE 


ADA 

-Rl 

5800 

ISOURCE 


STfl 

Resu 1 t 

5810 

ISOURCE 


LDfl 

-2 

5820 

ISOURCE 


LDB 

-147155B 

5830 

ISOURCE 


JSM 

Rel_math 

5840 

ISOURCE 


ISZ 

Count 

5850 

ISOURCE 


LDfl 

-3 

5860 

ISOURCE 


CPA 

Count 

5870 

ISOURCE 


JMP 

* + 2 

5880 

5390 

ISOURCE 

ISOURCE 

i 

JMP 

Get_f req 

5900 

ISOURCE 


LDfl 

-Array 1 

5910 

ISOURCE 


LDB 

-768 

5920 

ISOURCE 

Cont i nue : 

CLR 

16 

5930 

ISOURCE 


ADA 

-16 

5940 

ISOURCE 


DSZ 

B 

5^50 

5960 

ISOURCE 

ISOURCE 

i 

JMP 

Cont i nue 


loop contain* operations which ar* 


Get the first word of the column of 
the data array which contains the 
velocity component for which we want 
to get the range. 


Mask and rotate to get the four 
bits containing the range. 


Subtract from 15 to get the 
actual range. 

Store the actual range in RangeN 
so that we can transfer it to the 
main program 

Use the range to find the power 
of two needed for the divisor. 


Convert the power of two into a 
real number. 


Divide 3.2E4 by the appropriate 
power of two, using BCD math. 


1 Decide whether to put the result 
1 in Rl, R2, or R3, depending on Count 


• Now, finally, call the utility to 
1 perform the division. 

1 Increment and check Count so as 
1 to follow the loop three times. 


1 Zero out the entire count-to- 
1 velocity conversion table so that 
1 it must be recalculated for each 
1 point. (This must be done if the 
■ mixing frequency or ranges are 
' changed between counts.) 
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5970 

ISOURCE 

LDA 

■Su 

Set initial values of Su, Sv, 

5980 

I SOURCE 

LDB 

-15 

Su, Suu, etc. to zero. 

5990 

ISOURCE 

CUv: CLR 

4 


6000 

ISOURCE 

ADA 

-4 


6010 

ISOURCE 

DS2 

B 


6020 

ISOURCE 

jmp 

Clear 


6030 

ISOURCE 

1 



6040 

ISOURCE 

LDfl 

-Bin u 

Clear the areas of memory which 

6050 

ISOURCE 

LDB 

-192~ 

util be used to ho<d the number of 

6060 

ISOURCE 

Cl car 2: CLR 

16 

samples per count for the histograms. 

6070 

ISOURCE 

ADA 

-16 


6080 

ISOURCE 

DS2 

B 


6090 

ISOURCE 

JMP 

Cl ear 2 


6100 

ISOURCE 

i 



6110 

ISOURCE 

LDfl 

-Thet a 

Convert Theta from degrees to 

6120 

ISOURCE 

STfl 

Oper 1 

radians using the Rel_math 

6130 

ISOURCE 

LDfl 

-Rad~ 

utility. 

6140 

ISOURCE 

STfl 

0per_2 


6150 

ISOURCE 

LDfl 

-Xvar 


6160 

ISOURCE 

STfl 

Result 


6170 

ISOURCE 

• LDfl 

-2 


6180 

ISOURCE 

LDB 

-147155B 

. 

6190 

ISOURCE 

JSM 

Rel_math 


6200 

ISOURCE 

i 



6210 

ISOURCE 

LDfl 

-Xvar 

Find the sine and cosine of 

6220 

ISOURCE 

STfl 

Oper 1 

Theta, and store them in the 

6230 

ISOURCE 

LDfl 

-S i n~ 

locations Sin and Cos, respect i vel y. 

6240 

ISOURCE 

STfl 

Result 


6250 

ISOURCE 

LDfl 

-1 


6260 

ISOURCE 

LDB 

-34213B 


6270 

ISOURCE 

JSM 

Rel_math 


6280 

ISOURCE 

LDfl 

-Cos 


6290 

ISOURCE 

STfl 

Resul t 


6300 

ISOURCE 

LDfl 

-1 


6310 

ISOURCE 

LDB 

-34224B 


6320 

ISOURCE 

JSM 

Re l_math 


6330 

ISOURCE 

i 



6340 

ISOURCE 

i The Itjop defined by Get int calculates intermediate values 

6350 

ISOURCE 

' used in converting counts 

to velocities. These values are the 

6360 

ISOURCE 

• same -'or al 1 

samples in a 

point, so they can be calculated 

6370 

ISOURCE 

1 separately. The loop calculates <Nub-NumixN) and <Mill*DfN>, 

6380 

ISOURCE 

• where N is 1, 

2, and 3. 


6390 

ISOURCE 

LDfl 

-3 


6400 

ISOURCE 

STfl 

Count 


6410 

ISOURCE 

Get i nt : LDfl 

Count 


6420 

ISOURCE 

ADA 

— 1 


6430 

ISOURCE 

SHL 

2 


6440 

ISOURCE 

STfl 

Offset 


6450 

ISOURCE 

LDfl 

-Nub 

Find Numi xN-N jb-Num i xN. 

6460 

ISOURCE 

STfl 

Oper_l 


6470 

ISOURCE 

LDfl 

-Numi xl 


6480 

ISOURCE 

ADA 

Offset 


6490 

ISOURCE 

STfl 

0per_2 


6500 

ISOURCE 

STfl 

Resu 1 t 


6510 

ISOURCE 

LDfl 

-2 


6520 

ISOURCE 

LDB 

-146717B 


6530 

ISOURCE 

JSM 

Rel math 


65 40 

ISOURCE 

LDfl 

-Df 1 

Find DfN-Mi 1 1 vDfN. 

6530 

ISOURCE 

RDfl 

Offset 


6560 

ISOURCE 

STfl 

Oper_l 
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6570 

ISOURCE 


STfl 

Result 


£580 

ISOURCE 


LDfl 

-Mill 


£590 

ISOURCE 


STfl 

Oper 2 


££00 

ISOURCE 


LDfl 

-2 


££10 

ISOURCE 


LDB 

-147037B 


6£20 

ISOURCE 


JSM 

Rel_math 


££30 

ISOURCE 


DSZ 

Count 


££40 

ISOURCE 


JMP 

Get_i nt 


££S0 

ISOURCE 

f 




6££0 

ISOUPCE 

1 

This itction 

f 1 Iters the 

data to remove counts which are 

££70 

ISOURCE 

1 

txctiiivtly far from the 

mean. First the section creates 

££80 

ISOURCE 

I 

a lit*, of the 

number of 

times each particular count appears 

££90 

ISOURCE 

1 

The list consists of 3 1024 word arrays (called Bin_u, Bin_v, 

£700 

ISOURCE 

1 

and Bin w> where the address of each word corresponds to the 

£710 

ISOURCE 

1 

count, and the value of 

the word indica.es the numoer of times 

£720 

ISOURCE 

1 

that particular count has appeared. 

£730 

ISOURCE 

1 




£740 

ISOURCE 


LDfl 

-0 


£750 

ISOURCE 


STfl 

Count 


£7€0 

ISOURCE 


LDfl 

-5 

• Set Check to determine whether the 

£770 

ISOURCE 


STfl 

Check 

* count or 1024 minus the count is to 

£780 

ISOURCE 


LDfl 

Ns 

1 be used. 

£790 

ISOURCE 


STfl 

I 


£800 

ISOURCE 

1 




£810 

ISOURCE 

F i 

11 bin: LDfl 

Count 


£820 

ISOURCE 


LDB 

Ns 

1 Figure out which element of the 

£830 

ISOURCE 


MPY 


1 data array we want to pick up. 

£840 

ISOURCE 


ADR 

I 


£858 

ISOURCE 


ADA 

— 1 


£8£0 

ISOURCE 


STfl 

El ementd 


£370 

ISOURCE 


LDA 

- 1 -it 

' Get a raw datum from the data array 

£880 

ISOURCE 


LDB 

-flrrayd 


£890 

ISOURCE 


JSM 

Get element 


€900 

ISOURCE 


LDfl 

Int" 


£910 

ISOURCE 


LDB 

-1023 

1 Strip off the first six 

€920 

ISOURCE 


AND 

B 

1 bits of the raw data word. 

£930 

ISOURCE 


LDB 

Count 

1 See if Count = Check. 

£940 

ISOURCE 


TCB 



£950 

ISOURCE 


ADB 

Check 


£9€0 

ISOURCE 


S2B 

Sk i p 

• If true, use the modified data 

£970 

ISOURCE 


TCfl 


< word as an index. If not, use 

£930 

ISOURCE 


ADA 

-1024 

1 1024 minus the data word. 

£990 

ISOURCE 

Skip: STfl 

Int 

1 Store the count we have gotten. 

7000 

ISOURCE 


LDfl 

Count 


7010 

ISOURCE 


LDB 

-1024 


7020 

ISOURCE 


MPY 



7030 

ISOURCE 


ADA 

Int 

i Go to the address in the array 

7040 

ISOURCE 


ADA 

— 1 

• corresponding to the value of 

7050 

ISOURCE 


ADA 

-Bin u 

• the count. 

70£0 

ISOURCE 


132 

A, I 

1 Increment the word by one, 

7070 

ISOURCE 


1*52 

Count 

1 indicating that one more count 

7080 

ISOURCE 


LDfl 

—3 

i with that value has been read. 

7090 

ISOURCE 


fljfl 

Count 


7100 

ISOURCE 


SZfl 

*+2 


7110 

ISOURCE 


JUP 

Fi 1 1 _bi n 


7120 

ISOURCE 


S r fl 

Count 


7130 

ISOURCE 


DsZ 

I 


7140 

ISOURCE 


;r.p 

F i 1 1 _b i n 


7150 

ISOURCE 

1 




7 1 £0 

ISOURCE 

1 

When '5dev>0, 

use this section to filter the data. 
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•770 

I SOURCE 

1 



•780 

ISOURCE 


LDR 

-Ns 

’790 

ISOURCE 


STfl 

Oper_l 

•360 

ISOURCE 


LDR 

-Xvar 

’310 

ISOURCE 


STfl 

Resul t 

•320 

ISOURCE 


JSM 

Int_to_rel 

•830 

ISOURCE 

1 



’840 

ISOURCE 


LDB 

Count 

•830 

ISOURCE 


SBL 

2 

’860 

ISOURCE 


RDB 

-flvg_u 

•870 

ISOURCE 


STB 

Oper 1 

’880 

ISOURCE 


STB 

Result 

’890 

ISOURCE 


LDB 

-Xvar 

’900 

ISOURCE 


STB 

Oper 2 

’910 

ISOURCE 


LDfl 

-2 

’920 

ISOURCE 


LDB 

-147135B 

.’930 

ISOURCE 


jsh 

Rel_math 

.’940 

ISOURCE 

1 



7930 

ISOURCE 


LDfl 

Count 

7960 

ISOURCE 


ADA 

-1 

7970 

ISOURCE 


STfl 

Count 

7980 

ISOURCE 


ADA 

—3 

7990 

ISOURCE 


SZfl 

»+2 

3000 

ISf.URCE 


JMP 

Get_avg 

3010 

ISOURCE 

1 



3020 

ISLURCE 


LDfl 

-0 

3030 

ISOURCE 


STfl 

Count 

3040 

ISOURCC 


LDfl 

-Zero 

3036 

ISOURCE 


LDB 

Count 

8060 

ISOURCE 


SBL 

2 

8070 

ISOURCE 


STB 

Offset 

3080 

ISOURCE 


RDB 

-Dev u 

8090 

ISOURCE 


XFR 

4 

8100 

ISOURCE 


LDfl 

Offset 

3110 

ISOURCE 


ADA 

-flvg_u 

8120 

ISOURCE 


STfl 

Int 

8130 

ISOURCE 


LDfl 

Count 

8140 

ISOURCE 


ADA 

-1 

8130 

ISOURCE 


LDB 

-1024 

8160 

ISOURCE 


MPY 


8170 

ISOURCE 


ADA 

■8i n_u 

8180 

ISOURCE 


STfl 

Address 

8190 

ISOURCE 


LDfl 

-1023 

8200 

ISOURCE 


STfl 

I 

8210 

ISOURCE 

i 



8220 

ISOURCE 

count : 

LDfl 

-I 

8230 

ISOURCE 


STfl 

Oper_l 

8240 

ISOURCE 


LDfl 

-Xvar 

8250 

ISOURCE 


STfl 

Resul t 

8260 

ISOURCE 


JSM 

Int_t o_rel 

8270 

ISOURCE 

i 



8280 

ISOURCE 


LDfl 

-Xvar 

8290 

ISOURCE 


STfl 

Oper_l 

8300 

ISOURCE 


LOR 

Int 

8310 

ISOURCE 


STfl 

0per_2 

8320 

ISOURCE 


LDfl 

-Yvar 

8330 

ISOURCE 


STfl 

Resul t 

8340 

ISOURCE 


LDfl 

-2 

8330 

ISOURCE 


LDB 

-146717B 

8360 

ISOURCE 


JSM 

Pel math 


' Convert the number o * samples to 
i a real number. 


i Divide the sum stored in flvg_u by 
! the number of samples to get the 
I average value of the counts. 


' Get the standard deviation of each 
' set of counts using the average and 
• the information in the Bin_u arrays. 

' First, set the standard deviation -0. 


i Get a count and convert it to a 
1 real number. 


' Subtract the average count from the 
i count we just got. 
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8370 

1 SOURCE 1 



8330 

ISOURCE 

LDA 

■ Yvar 

8390 

ISOURCE 

STA 

Oper_l 

8400 

ISOURCE 

STA 

0p*r_2 

8410 

ISOURCE 

LDA 

■ Xvar 

8420 

ISOURCE 

STA 

Resul t 

8430 

ISOURCE 

LDA 

-2 

8440 

ISOURCE 

LDB 

-147037B 

84S0 

ISOURCE 

JSM 

Rel_r»ath 

8460 

ISOURCE ' 



8470 

ISCURCE 

LDA 

Address 

8480 

ISOURCE 

STA 

Oper_l 

8490 

ISOURCE 

LDA 

■Yvar 

8S00 

ISOURCE 

STA 

Resul t 

8510 

ISOURCE 

JSM 

Int_t o_re 

8520 

ISOURCE ' 



8530 

ISOURCE 

LDA 

■Xvar 

8540 

ISOURCE 

STA 

Oper_l 

8550 

ISOURCE 

LDA 

■Yvar 

8560 

ISOURCE 

STA 

0per_2 

8570 

ISOURCE 

STA 

Result 

8580 

ISOURCE 

LDA 

-2 

8598 

ISOURCE 

LDB 

■147037B 

8600 

ISOURCE 

JSM 

Rel_math 

8610 

ISOURCE • 



8620 

ISOURCE 

LDA 

■Yvar 

8630 

ISOURCE 

STA 

Oper_l 

8640 

ISOURCE 

LDA 

■ Dev u 

8650 

ISOURCE 

ADA 

Offset 

8660 

ISOURCE 

STA 

0per_2 

8670 

ISOURCE 

STA 

Resul t 

8680 

ISOURCE 

LDA 

■2 

8690 

ISOURCE 

LDB 

■ 1 4672 1 B 

8700 

ISOURCE 

JSM 

Re1_ math 

8710 

ISOURCE • 



8720 

ISOURCE 

DSZ 

Address 

8730 

ISOURCE 

DSZ 

I 

8740 

ISOURCE 

JM° 

Dev_count 

8750 

ISOURCE ' 



8760 

ISOURCE 

LDA 

■Ns 

8770 

ISOURCE 

STA 

Oper_l 

8780 

ISOURCE 

LDA 

■Xvar 

8790 

ISOURCE 

STA 

Result 

8800 

ISOURCE 

JSM 

Int_to_re 

8810 

ISOURCE ' 



8820 

ISOURCE 

LDB 

Offset 

8830 

ISOURCE 

ADB 

■Dev_u 

8840 

ISOURCE 

SIB 

Oper_l 

8850 

ISOURCE 

STB 

Resul t 

8860 

ISOURCE 

LDB 

■ Xvar 

8870 

ISOURCE 

STB 

Oper 2 

8880 

ISOURCE 

LDA 

-2 

8890 

ISOURCE 

LDB 

■147155B 

8900 

ISOURCE 

JSM 

Rel_math 

8910 

ISOURCE ' 



8920 

ISOURCE 

LDA 

■ Dev u 

8930 

ISOURCE 

ADA 

Offset 

8940 

ISOURCE 

STA 

Oper_l 

8950 

ISOURCE 

STA 

Resu 1 t 

8960 

ISOURCE 

LDA 

■1 


1 Squire the difference between the 
1 tvertgt and the count we got to 
' get the deviation from the mean. 


• Use the Bin_u arrays to find out 

' how many counts there are with the 
■ value ue have chosen, and convert 

* that number to a real number. 


• Multiply the number of counts by 
■ the deviation from the mean. 


Add the product to Dev_u, which 
currently contains a running sum of 
the number of counts times the 
deviation from the mean of each 
count . 


• Convert the number of samples to 
' a real number. 


! Divide the sura stored in Dev_u by 

• the number of samples to get the 

• variance of the counts. 


I Take the square root of the 
• variance to get the standard 
' deviation. 
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8970 

ISOURCE 


LDB 

•31450B 

8980 

ISOURCE 


JSM 

R*l_math 
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ISOURCE 
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9000 

ISOURCE 


LDfl 

Count 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 


LDA 
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ISOURCE 
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ISOURCE 

f 1 1 “ r : 
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ISOURCE 


ADA 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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Address 

9150 

ISOURCE 


LDA 

-1023 

9160 

ISOURCE 


STA 

I 

9170 

ISOURCE 
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ISOURCE 


LDA 

Count 

9190 

ISOURCE 


SAL 
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9200 

ISOURCE 


ADA 

-Dev_u 

9210 

ISOURCE 


STA 

Op*r_l 
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ISOURCE 


LDA 

-Sdev 
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ISOURCE 


STfl 

0p*r_2 

9240 

ISOURCE 


LDA 

-Xvar 

9250 

ISOURCE 


STA 

Resul t 
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ISOURCE 
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-2 

9270 

ISOURCE 


LDB 

•1470378 

9280 

ISOURCE 


JSM 

R*l_math 
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ISOURCE 
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ISOURCE 


LDfl 

-Xvar 

9310 

ISOURCE 


STfl 

Oper_l 
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ISOURCE 


LDfl 

-Int 
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ISOURCE 
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Result 
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ISOURCE 
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ISOURCE 
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ISOURCE 


LDfl 

Count 
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ISOURCE 
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2 

9380 

ISOURCE 
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9390 

ISOURCE 
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Op*r_l 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 


LDB 

Offset 
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ISOURCE 
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ISOURCE 
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I 
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ISOURCE 
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B 
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ISOURCE 
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*+2 
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ISOURCE 


TCfl 
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ISOURCE 


ADA 
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ISOURCE 


SAP 

*+3 
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ISOURCE 


LDfl 

Address 

9530 

ISOURCE 


CLR 

1 

9540 

ISOURCE 


DSZ 

Address 

9550 

ISOURCE 


DSZ 

I 

9560 

ISOURCE 


JMP 

Fi 1 tr dev 


• Now us* the standard deviation to 
1 filter out all the counts whose 
1 value is more than Sdev standard 
1 deviations away from the mean. 
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ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 
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ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 
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ISOURCE 
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ISOURCE 


ISZ Count 
LOR --3 
ADR Count 
SZR *+2 
JMP Dev f i 1 tr 


1 Co through this whole process 
! of checking the deviation an d 
* filtering the data three times, 
' once for eac.i channel. 


The loop Begin is performed three times, once for U, V, 
and W. Each time througn the data array Is read, a count (a raw 
datum) is taken from it and converted into a velocity. The vel- 
ocity is stored in U, V, or 14 depending on whether this is the 
first, second, or third iteration of the loop. The first time a 
particular count is encountered, the velocity correspond! ng to 
it is calculated using the intermediate values found in Get_*nt, 
and the velocity is stored in a table. If that count is found 
again in the data array, the correspond) ng velocity is looked 
up rat her than being calculated again. 


Set begin: LDA 
STR 

' LOR 

STR 
STR 

i 

Begin: LOR 


Strai ght : 


>0 

Count 

Ns 

1 

Ngs 

Count 

B 

2 

Offset 

Ns 


--1 

El em*ntd 

-Irit 

-Rrrayd 

Get_e 1 erent 

I nt 

-1023 

B 

Count 

Check 
Strai ght 

*1024 

Int 

Count 

-1024 

Int 
— 1 

-Bi n_u 
fl, I 

Goodcount 

Court 

Ngs 

I 

Beg i n 


1 Figure out which element of the 
1 data array we want to pick up. 


• Get a raw datum from the data array. 


1 Strip off the first six 
1 bits of the raw data word. 
1 See if Count - Check. 


1 If true, use the modified data 
i word as an i idex. If not, use 
' 1024 minus the data word. 

• Store the count we have gotten. 

1 Now look in the appropriate part 
1 of the arrays created by the filter 
1 section. 


If the word correspondi ng to the 
count is zero, the count was 
filtered out and should be ignored. 
Ignore these three counts, go back 
and do the ne <t three. 
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10170 

ISOURCE 

Goodcount : 

LOR 

Count 

10180 

I SOURCE 


LOB 

-4096 

10190 

ISOURCE 


MPY 


10200 

ISOURCE 


LOB 

Int 

10210 

ISOURCE 


ROB 

— 1 

10220 

ISOURCE 


SBL 

2 

10230 

ISOURCE 


ROR 

B 

10240 

ISOURCE 
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-1 
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ISOURCE 


ADR 

-flrrayl 
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ISOURCE 


STfl 

Address 
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ISOURCE 


LOB 
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ISOURCE 


XFR 
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10290 

ISOURCE 


LOR 

Int 2 
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ISOURCE 


S2R 

Calcui ate 
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ISOURCE 


JMP 

Over 

10320 

ISOURCE 
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LOR 
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ISOURCE 
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Oper_l 

10340 

ISOURCE 


LOR 

-Yvar 

10350 

ISOURCE 

• 

STR 

Result 
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ISOURCE 


JSM 

Int_t o_rel 
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ISOURCE 


STB 

Oper 2 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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-147155B 

104S0 

ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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10710 

ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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Oper 1 
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ISOURCE 


LOR 

-Of l - 
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ISOURCE 
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Off*-* 
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• Now, use Count to find out 
' which lookup table array 

1 we want to use, and us* the 
1 count we got from the data array 

• to find exactly where in the 
1 table we want to go. 


• If that table entry is zero, 
< calculate a velocity for it. 


• Convert the count into 

• a real number. 


1 Tivide the range we found 
1 earlier by the count to get a 
1 frequency. 


i Find <Nub-Numi xN)-FrequencyN. 


1 If we are calculating U, reverse 
• the sign of < <Nub-Numi xN)-FrequencyK 
1 so as to reverse the sign of U. 

1 Leave V and Wv alone. 


i Find Vel oc i ty-< (Nub-Numi xN> 

• -Frequ*ncyN>*(Mi 1 1 *DfN> 

1 and store in a place m the 

• lookup table corresponding to 


i'f t 


/ 


/ 


t 
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10770 

ISOURCE 


STfl 

Qper_2 1 

the data count . 

10730 

ISOURCE 


LDfl 

Address 


10790 

ISOURCE 
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10800 

ISOURCE 


STfl 

Result 


10810 

ISOURCE 


LDfl 

-2 


10820 

ISOURCE 


LDB 

-147037B 


10830 

ISOURCE 


JSM 

Rel_m ath 
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ISOURCE 
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10838 

ISOURCE 

Ovtr! LDfl 

Address 1 

Transfer the velocity from 

10860 

ISOURCE 


ADA 

— 1 ! 

the lookup table to U, V, or 

10870 

ISOURCE 


LDB 

«U 1 

as appropriate. 
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ISOURCE 


ADB 

Offset 


10890 

I oOURCE 


XFR 

4 


10900 

ISOURCE 
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10910 

ISOURCE 


ISZ 

Count 


10920 

ISOURCE 


LDfl 

■ -3 ' 

Have 0, V, and 14 al 1 

10930 

ISOURCE 


ADA 

Count 1 

been calculated* 7 

10940 

ISOURCE 


SZfl 

♦+2 1 

If not, go back again. 
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ISOURCE 


JMP 

Begin 


18960 

ISOURCE 


STfl 

Count 1 

If so, set Count ■ 0. 
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ISOURCE 
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10980 

ISOURCE 

1 

Now u# convert the 14 we have obtained (which is measure 

10990 

ISOURCE 

1 

angle Theta to the V-axis> 

to the 14 we want (which shou 

11000 

ISOURCE 

1 

be me asured 

at an angle of 

90 degrees to the V-axis>. 

11010 

ISOURCE 

1 

Thus -'ind 14- 

(l4v-V*Cos (ThetaJJ-'Si n< Thet a) . 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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I SOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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Resul t 
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ISOURCE 


STfl 

0per_2 
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ISOURCE 
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ISOURCE 
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Oper_l 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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Oper_l 
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ISOURCE 
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ISOURCE 
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I SOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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Now tike running sums of U, 

V, 14, and several products 

11320 

ISOURCE 

1 

of th*s# velocities. The sums are taken using the ut i li 

11330 

ISOURCE 

1 

“Add". The sums are calculated in an unusual sequence i 

11340 

ISOURCE 

1 

order to reduce the number 

of program steps needed to c 

11350 

ISOURCE 

1 

culat* them. 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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-Suu 

' Find 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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1 Find 
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ISOURCE 
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Oper_i 



ISOURCE 
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ISOURCE 
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Result 



ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 

LOR 
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ISOURCE 
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Resul t 
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ISOURCE 

LDfl 
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ISOURCE 
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ISOURCE 
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Rel_m*th 
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ISOURCE ' 
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ISOURCE 

LDfl 

•Suvv 

' Find 
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ISOURCE 
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Resul t 
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ISOURCE 
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ISOURCE ' 
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ISOURCE 
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■ Find 
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ISOURCE 
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Open 1 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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Resul t 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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Rel_m*th 
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ISOURCE ' 





12170 

ISOURCE 

LDfl 
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1 Find 

Suww*Suuw+<U*V*W) 

12180 

ISOURCE 

LDB 

Resul t 



12190 

ISOURCE 

JSM 

Add 
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ISOURCE ' 
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ISOURCE 

LDfl 
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! Find 

V#V 
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ISOURCE 
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Oper_l 
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ISOURCE 
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0per~2 
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ISOURCE 
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ISOURCE 
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Result 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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Rel_m*th 
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ISOURCE ' 





12300 

ISOURCE 

LDfl 

*Suu 

1 Find 

Svv*Suv+< V*V> 

12310 

ISOURCE 

LDB 

Result 
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ISOURCE 

JSM 

Add 



12330 

ISOURCE i 
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12340 

ISOURCE 

LDfl 
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' Find 

Sw*Sw+W 
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ISOURCE 

LDB 
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ISOURCE 

JSM 

Add 
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ISOURCE ' 
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ISOURCE 

LDfl 
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' Find 

u*w 
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ISOURCE 
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Oper_l 
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ISOURCE 

LDfl 

*Xv*r 
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ISOURCE 
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Resul t 



12420 

ISOURCE 
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-2 
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ISOURCE 
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ISOURCE 
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ISOURCE ' 
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ISOURCE 
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' Find 
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ISOURCE 
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Resul t 
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ISOURCE 
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Add 
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ISOURCE ' 
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ISOURCE 
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' Find 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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Oper_l 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 


LDA 

■2 
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ISOURCE 
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ISOURCE 
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Rel_math 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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12750 

ISOURCE 
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ISOURCE 
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ISOURCE 
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Oper_l 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 


JSM 

Add 
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ISOURCE 
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ISOURCE 
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1 Cont 
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ISOURCE 
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1 sums 

12910 

ISOURCE 
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12920 

ISOURCE 

• Now 

place the finished 

summs in 

12930 

ISOURCE 

1 the 

BASIC program has . 

access to 

12940 

ISOURCE 

1 the 

BASIC pr% 

gram. 
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ISOURCE 
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12960 

ISOURCE 
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ISOURCE 


JSM 

Put_value 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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ISOURCE 
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U*U 


Sww»Sww+(M*W> 


u*w*w 


Suuw=»Suuw + <U*W*W> 


nue te calculate running 
until out of samples. 

the COMMON region so that 
them, and then return *o 


— v 


4 


13170 

ISOURCE 


LOB 

■Suu_par 

13180 

ISOURCE 


JSM 

Put_ua1 ue 

13190 

ISOURCE 


LOR 

»Suu 

13200 

ISOURCE 


LOB 

■Suu_par 

13210 

ISOURCE 


JSM 

Put_ua1 ue 

13220 

ISOURCE 


LOR 

■ Suu 

13230 

ISOURCE 


LOB 

■Suu_par 

13240 

I'. URCE 


JSM 

Put_ua1 ue 

132S0 

ISOURCE 


LOR 

■Suw 

13260 

ISOURCE 


LOB 

■Suu_par 

13270 

ISOURCE 


JSM 

Put_ual ue 

13280 

ISOURCE 


LOR 

■Suu 

13290 

ISOURCE 


LOB 

■Suu_p ar 

13300 

ISOURCE 


JSM 

Put_val ue 

13310 

ISOURCE 


LOR 

■Suuu 

13320 

ISOURCE 


LOB 

■Suuu_par 

13330 

ISOURCE 


JSM 

Put_val ue 

13340 

ISOURCE 


LOR 

■ Suuu 

13350 

ISOURCE 


LOB 

■Suuu_par 

13360 

ISOURCE 


JSM 

Put_ual ue 

13370 

ISOURCE 


LOR 

■ Suuu 

43380 

ISOURCE 


LOB 

■Suuw_par 

13390 

ISOURCE 


JSM 

Put_ual ue 

13400 

ISOURCE 


LOR 

>Suuw 

13410 

ISOURCE 


LOB 

■Suwu_par 

13420 

ISOURCE 


JSM 

Put_ua> ue 

13430 

ISOURCE 


LOR 

■Suuu 

13440 

ISOURCE 


LOB 

■Suvu_par 

13450 

ISOURCE 


JSM 

Put_ual ue 

13460 

ISOURCE 


LOR 

■Ngs 

13470 

ISOURCE 


LOB 

■Ng*_par 

13480 

ISOURCE 


JSM 

Put_ual ue 

13490 

ISOURCE 

1 



13500 

ISOURCE 


LOR 

■ B» n u 

13510 

ISOURCE 


LOB 

■ 192~ 

13520 

ISOURCE 


STB 

I 

13530 

ISOURCE 


LOB 

■16 

13540 

ISOURCE 


STB 

Count 

13550 

ISOURCE 


LOB 

Cbin par 

13560 

ISOURCE 

Transfer: 

XFR 

16 

13570 

ISOURCE 


ADR 

Count 

13580 

ISOURCE 


ROB 

Count 

13590 

ISOURCE 


OSZ 

I 

13600 

ISO'JRCE 


JMP 

Transf er 

13618 

ISOURCE 


RET 

1 

13620 

ISOURCE 

i 



13630 

ISOURCE 


LIT 

200 

13640 

ISOURCE 

J 



13650 

ISOURCE 

1 The u% i 

1 1 t y> ' 

"fidd" is u 

13660 

ISOURCE 

• 



13670 

ISOURCE 

Add: 

IS7 

Ut t count 

13680 

ISOURCE 


STfl 

Oper_l 

13690 

ISOURCE 


STB 

0per_2 

13700 

ISOURCE 


STfl 

Resu 1 1 

13710 

ISOURCE 


LOR 

■2 

13720 

ISOURCE 


LOB 

■ 1 4672 1 B 

13730 

ISOURCE 


JSM 

Re 1 _mat h 

13740 

ISOURCE 


OSZ 

Ut 1 c ownt 

13750 

ISOUPCE 


RET 

1 

13760 

ISOURCE 


JSM 

Ut 1 end 


Transfer the contents of Bin_n 
to the common area without using 
the siou HP-supplied external 
subrout 1 nes . 


to add up the running sums. 
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TABLE 1 

MACRODYNE RANGE SETTINGS 



RANGE 


0 

32 GHz 

64 MHz 

31.28 MHz 

1 

16 ” 

32 " 

15.64 

9 f 

2 

8 " 

16 " 

7.20 

V 9 

3 

4 ” 

8 " 

3.91 

9 9 

4 

2 M 

4 " 

1.95 

9 9 

5 

1 ” 

2 •* 

977 kHz 

6 

500 MHz 

1 M 

488 

9 9 

7 

250 " 

500 kHz 

244 

9 9 

8 

125 " 

250 " 

122 

9 9 

9 

62.5 " 

125 ” 

61 

9 1 

10 

31.25 " 

62.5 M 

30 

9 9 

11 

15.62 ” 

31.25 ” 

15 

9 9 

12 

7.81 " 

15.62 " 

7.6 

9 9 

13 

3.90 ” 

7.81 ,r 

3.8 

9 9 

14 

1.95 " 

3.90 ” 

1.9 

9 9 
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c tULUti ruud 
BEAM MATRIX 
TO TRANSMITTING 
TABLE 

ROTATING 

PRISM 

BEAM 

DISPLACER 

BEAM 

STEERING 

MODULES 


Q 


COLLIMATOR 


BRAGG CELL 
(BLUE) 

(14) 

BRAGG CELL 
(GREEN) 

(13) 

BEAM SPLITTER 
(BLUE) 

(12) 

POLARIZATION 

ROTATOR 

(BLUE) 

(11) 

BEAM DISPLACER 
(BLUE) 

(10) 

BEAMSPLITTER 

(GREEN) 

(9) 

POLARIZATION 

ROTATOR 

(GREEN) 

(8) 


GREEN 514 3 nm4 ♦ BLUE 488 nm 


MIRROR 

a 


ATTENUATOR 

(3) 


DISPERSION 
PRISM 
(4) 


MIRROR 

(7) 


MIRROR 

(61 


COLOR SEPARATOR BOX 


Schematic of the optics table layout 








Fig. 2 Schematic cf the transmitting optics. 




















**4 

CJ 



Fig. 4 


Schematic of the receiving optics 












































RESETEVENT 


EVENT PULSE 


EVENT 

I SYNCHRONIZER T j M e/stATUS > ,DATA 

L__ — r — 


DIGITAL 
DATA ✓ 

'16 TO 

HP 9845 
r CPU 

handshake 


A/D 

digitized 

CONVERTER 

ANALOG ? 


Fig. 7 


Computer interface simplified block diagram. 





i 




Fig. 8 


Block diagram for data acquisition program. 
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Fig. 9 


Block diagram for data reduction program. 
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SYNC PULSES FROM COUNTER 


DIGITAL MULTIPLEXER INPUTS 16 BIT WORDS 







’■NTtHfACe 


Fig. 10 NASA LDV-A/D computer interface connections 

and settings. 
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I 


HONEYCOMB 






Fig. 11 Experimental rig. 
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2 


(b) LDV neasurenants 


Fig. 12 


Secondary valocity plot#. 
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